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ABSTRACT 


This thesis is a continuation of previous work to apply 
modern multiparameter estimation techniques to the problem 
of estimating attrition rates for a large number of small 
inventory cells in manpower planning models used by the U.S. 
Marine Corps. The main advances involve the promising 
introduction of empirical Bayes (non-constant shrinkage) 
techniques, recognition of the non symmetric nature of the 
errors with a response to this, and some insight into all 
aggregation plans that should help provide greater stability 
for the estimation methods. In addition, the roles of some 


middle level methodological choices are explored. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs 
developed in this research may not have been exercised for 
all cases of interest. while every effort has been made, 
within the time available, to ensure that the programs are 
free of computational and logic errors, they cannot be 
considered validated. Any application of these programs 


without additional verification is at the risk of the user. 
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A. PURPOSE 

MPs pDOobe or project OPUS (Orfticer Planning and Utility 
System), Headquarters, U.S. Marine Corps, in 1985, requested 
assistance in exploring new methods of generating manpower 
loss rates in order to improve upon the currently used 
method of force rate analysis. In response, students and 
faculty in the Operations Research department at the Naval 
Postgraduate School began exploring the use of modern 
multiparameter estimation techniques for τηλ Ес Ти. 
Special emphasis is placed on the small celi problem, i.e., 
categories of officer skill, grade and length of service 
which have low inventory figures. ხა ი 1C3) race 
estimators for small cells are unstable and a large number 
of these cells exist. 

MTS paper builds upon previous work on this problem. 
The main advances are to append new measures of 
effectiveness as requested by the Navy Personnel Research 
and Development Center  (NPRDC), introduce a class of 
empirical Bayes estimators, and explore the effects of some 
middle level choices in applying the new and existing 
methods. The ultimate goal 15 to refine the techniques 
presented here in order to validate a clear policy for 


predicting loss rates. 


B. BACKGROUND 

For an introduction into Marine Corps policy concern 
manpower planning, the reader is referred to a thesis 
submitted in September 1985, at the Naval Postgraduate 
School, by "Major D.D δι er τε 11, who presents a 
detailed background into the Marine Corps officer structure 
and the manpower planning process. Tucker also provides 
basic attrition rate theory and calculates overall attrition 
rates in several different formats as they pertain to the 
Marine Corps. As the aggregation rate begins to grow 
smaller with refinement, Tucker illustrates the irregular 
behavior of losses due to voluntary attrition and the small 
numbers of losses obtained in such aggregates. Tucker [Ref. 
1:p. 50] introduces the James-Stein technique of loss rate 
estimation as a shrinking of individual cell averages toward 
a grand average in order to reduce the risk, or squared 
differences of forecast and actual values, the goal being an 
improvement over the classical maximum likelihood estimator, 
referred to from now on as the MLE. Tucker chose the ranks 
of First Lieutenant and Lieutenant Colonel for evaluation 
using military occupational specialty (MOS) groups of combat 
Support and ground combat (three MOS's were selected for 
each group). These aggregation schemes were carried through 
other studies and are included for historical purposes in 


this paper. 


Another thesis submitted at the Naval Postgraduate 
от zn March 1986, by Major J.R. Robinson (Ref. 2] 
focused on a technique called limited translation of the 
James-Stein estimator. This particular model attempts to 
minimize the risk mentioned above by reducing the shrinkage 
of rate estimates for cells which fall outside a certain 
range of values centered on the grand mean. It was thought 
that shrinking all cells by the same rate toward this grand 
mean may be unwise since the attrition rates for those 
Barcjıcular cells, those farthest from the grand mean, may 
contain attributes which are different than the majority orc 
tne cells in the aggregate. Additionally, Robinson's work 
carries forward the study done by Tucker concerning MLE and 
James-Stein estimators. 

These feasibility studies were limited in scope due to 
the format of the summary date tape available (acquired from 
NPRDC), which allows only coarse cell definition. A cell 1s 
defined to be a cross-classification of the forty military 
occupational fields (OF), thirty-one lengths of service 
Mees), and ten grades for a total of 12,400 categories for 
manpower planning purposes [Ref. 2B. 1012 The data does 
mee. distinguish between limited duty or Cers, those 
officers specifically designated for limited duty within 
certain MOS's, and unrestricted officers (regular and 
reserve) [Ref. 1:р. 20]. These factors dictated a broad 


aggregation scheme that created stable cell inventories and 


allowed the common variance assumption crucial to their 
work. But some of the application models require more 
refined information along those lines and experiments with 
finer aggregation levels, necessary in some real cases, 
caused this assumption to be violated. 

Tucker and Robinson did conclude that the current schene 
of attrition rate prediction could be improved upon by their 
methods, although no dominant scheme was uncovered. Tucker 
and Robinson did identify problem areas in their studies 
that required additional work. Tucker [Ref. 1:p. 71] cited 
a need for a better aggregation method to produce a more 


4 


homogeneous attrition behavior in all cells. Also, small 
probabilities of loss within a cell were not dealt with 
successfully. Robinson [Ref. 2:0. 32] revealed the small 
cell problem and the inability to normalize the cell means 
or stabilize the variance with a data transformation. The 
recommendetions of the previous studies highlighted the need 
for a more refined data tape with current information to 
include, among others, full Military Occupational Specialty 
(MOS) information, grade separation to include 


regular/reserve status, promotion zone data and breakdown of 


attritions by type. 


er PROGRESS 
The prcject of loss rate estimation, sponsored by NPRDC, 
is currently moving into the implementation stage. An 


operational data tape is now available which includes 


4 


detailed information on Marine Corps officers for the years 
1977-86. Work has begun at the Naval Postgraduate School to 
break out the data in a useful format and its availability 
in now nearing completion. 

In addition to the results stated earlier, Tucker and 
Robinson provided loss rate estimates for small cells with 
no attrition, i.e., MLE equal to zero, and risk effects of 
the different schemes studied. These results provided the 
springboard for this study and the requirements for 
additional future studies. 

Captain R.W. Larsen [Ref. 3] has developed a promising 
aggregation scheme based on cluster analysis, 1.е., а 
classification scheme to aggregate cells which reflect a 
greater degree of homogeneity in attrition rates than are 

t 
allowed using current aggregation methods. Table 1 displays 


the resulting scheme, given in year of current service (YCS) 


groups proposed for implementation by Larsen. 


Table 1. AGGREGATION METHOD PROPOSED BY LARSEN 





MOS Category Bounded YCS Groups 
Fixed-Wing Pilots (1-6,8-19) (7) (20-25) (26) 
Rotary- Wing Pilots (1-5,8-19) (6,7) (20-25) (26) 


Naval Flight Officers (1-5,8-19) (6,7) (20-25) (26) 
Lawyers (1-6,8-19) (7) (20-25) (26) 
All Else ~ (1-3,6-19) (4,5) (20-25) (26) 


The parentheses encompass YCS groups which behave similarly. 
In each case the grade is fized, i.e., one should aggregate 
over YCS before aggregating over grade. 

This information will allow for more successful 
application of the loss rate estimation techniques of 


interest here when the new data tape is provided. 


D. GOALS 

In an attempt to refine the methods of attrition rate 
estimation presented in the previous pilot  stud:es, the 
following goals were set for this paper: 


1. Modify the existing пе ου σσ ος 

a) use the refined data format, 

b) extend the work to the case of unequal or non 
constant variance in the small cell inventory 
problem, 

c) study alternate transformation inversion 
techniques, 

d) introduce and evaluate additional measures of 
effectiveness suggested by NPRDC to include, but 
not limited to, cross validation underages and 
overages of the deviation between forecasts and 
actuals, 

e) examine the graphical effect of shrinkage in the 
original scale. 


2. Introduce empirical Bayes estimation method in several 
forms for attrition rates for consideration as another 
option to solve the problem. 


οι ZZ ERBEN STONSTROFESTIMATTON METHODS 


ES GENERAL 


As stated in Chapter I, Tucker and Robinson published 


prot studies introdücing the James-Stein, limited 
translation James-Stein, Maximum Likelihood and Transform 
Scale Cell Average estimators as alternatives og the 


aggregate methods previously used for manpower planning in 
the Marine Corps. Their performance was promising, and this 
chapter will investigate ways to enhance their performance 
and possibly develop a dominant estimator por | future 
ir) tlon prediction. Additionally, NPRDC has expressed 
interest in additional measures of effectiveness for these 
estimators which may lead to a sharper direction for the 
qoal of producing a viable estimation policy. This chapter 
concentrates on the following issues concerning previous 


studies. 


B. NON CONSTANT VARIANCE 

For a cell having an inventory n, and attrition rate XL, 
the number of attritions y is described by the Binomial 
(n,p) distribution. The variance of the estimator of 
р = y/n is given by the familiar formula 


var (p) „Bit pl | 


bo 
pa 


The Freeman-Tukey double arcsine transformation was used 
by Tucker [Ref. 1:p. 55] and Robinson [Ret. 3:55. MN 
order to map the raw losses to the transformed scale. е 


formula is Given БУ 


თ 


x =(0.5 (n+ ο ο το ааа 1] (2... 
Те 2(y + 1) _ 
+ sine ci ο σαν К 


Note that n=n(t) may change with time and thus, so would x. 
The question of averaging over time before or after applying 
the Freeman-Tukey transformation is not discussed in this 
paper. One of the goals of using this transformation is to 
stabilize the variance of this estimate. In thE sS 
Tucker and Robinson, the assumption of normally distributed 
random variables with common variance was the setting in 
which the James-Stein estimator was expected to perform 
well. Indeed, for n moderate or large and p not too 
extreme, the variance of the Freeman-Tukey transformation is 
approximately one. But as discovered by Robinson, this 
assumption does not hold well when applied to cells with 
small inventory. A more careful look at the variance of the 
transformed data is taken here, with particular interest in 
the region of unstable variance where n and p are small. 
Some exploratory graphical work has revealed that the 


natural log of the variance of x, the transformed value, 


behaves quite well as a linear function of u( =E[X]) and p- 
1 in this unstable топе. This leads to the interpolatory 


Formula: 


ყა ე)” 568355 2953300 - 1). 999. (2. 3) 


ΙΙ 2.2 


oper limit, U - 2.2, 15 the value for which var(x) = 1. 
Since the least squares fitting process was applied to the 
Mear thms, u can not be allowed to fall to one or below. 
„“ Cლ67CXL limit, n = 1.001, is an arbitrary value that meets 
mers constraint. The coefficients were validated bv a 
ear regression on the natural log of this formula. The 
emma) Variance function for n= 1,2,3,4,5,7,10 appear in 
Figure 1, andthe fit is remarkably good for n > 3. The 
upper graph shows the region of special interest for this 
mper and this portion of tHe curve is fitted by the formula 
ШЕЕсаласіосп 2.3 ароте. The lower plot exhibits the outer 


tale effect yet to be explored. The small values of n(<3) 


S5 hot maintain the stable variance region, 1.е., var[x]=1, 
for long before falling back. As the value for n increases, 
the more stable the variance becomes. Table 2 is a partial 


display of the residuals for selected n values, fifteen 
equally spaced values in the range 1.001 = п s 2.2, computed 
by the linear regression formula above. For values of nm > 2 
the difference between actual and fitted values are 
considered acceptable for this work. 
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Table 2. 








RESIDUAL LISTING FOR SELECTED VALUES OF N 





N= 1 B N=4 N= 5 N= 10 
0.07191 ОО? 9.02155 0.017 0.008313 
0.08466 0.02569 0.002853 -0.001695 -0.01106 
0.1033 0.0239 -0.005448 -0.01049 -0.0206 
00122 0.03002 -0.00655 -0.01185 -0.02257 
0.1585 0.04068 -0.004026 -0.009566 -0.02067 
DELS 0.05353 0.0000676 -0.005727 -0.0166 
01552 0.06853 0.004027 -0.002103 -0.01291 
021527 0.05234 0.0074-7 0.0009541 -0.009981 
0.1414 0.0948] 0.00948 0.002557 -0.008299 
0.1209 0.1055 0.009825 0.002448 -0.008128 
0.08958 0.1129 0.008569 0.0004583 -0.009929 
0.04832 0.1169 0.004997 -0.005511 -0.01527 
-0.005254 OCDE 0.0001116 -0.00595 -0.0182 
-0.06914 0.113 -0.006798 -0.01654 -0.02463 
-0.1466 0.104 -0.02585 -0.03366 


-0.01534 


The boundaries for u, in the unstable region, as a 


function of n and p are provided in Figure 2. The formula 


derived was quite adequate for this study. 


C. ALTERNATE TRANSFORM INVERSION METHODS 

There are three issues in the inversion of the data from 
the transformed space. The first is that of selecting 
values of n, i.e., an average of the n(t)'s appearing in the 
is discussed later 


estimation year inventories. This issue 


in the chapter, with the development of a simulation for 


comparing mean inventory values. The second is that of the 


inversion of the sum of the two arcsine functions. The 
third involves the case of time averaging in the transformed 


scale and the question of what matching form of average 
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maventory should be used in the inversion. The third issue 
is not covered in this paper, and is still pending. 
Treatment of the second issue will be provided first. 

[re method used in the pilot studies of Tucker and 
КОШ соп, and the most transparent transform inversion 


formula is: 


д 9 A 
O 7 511 (д + 1/27) (264) 


Mime will be referred to as the Basic Inversion. The 


ООЛО а offered by Carter and Rolph [Ref. 4] is given by: 


| 
(4 (n/k) + 2) 


I) 
Ul 


δι = sin? (i) + | ET БЕ 
where Өң is х/(п + 1/2)!/? and 8, is the estimated amount of 
shrinkage in the transformed scale depending on the Bayes 
wedel used (constant or proportional prior), k is the number 
ells and n is inventory. 

БЕсатепсіу, Carter and Rolph invert prior to shrinking 
and the above formula provides a shrinkage in the original 
scale. Rewriting the above, it appears to shrink towards 


me 2/2, and is given by: 


a (2.6) 


) * віп2 (91) (1 წ... = 


E 28, 
Рі = 1/27) + 2 


An exact  Freeman-Tukey inversion formula developed by 
Miller [Ref. 5] was recently uncovered. It is exact in the 


sense that for a single cell, the formula given by: 


m3 


S(t) = 0.51 а. 


sin S9 στι Е уг улуг | (2.7) 
n 
where ს = x/(n + 1/2)!/? + n/2 and a is the average central 
inventory, returns the empirica Tate y/n. This will be 


referred to as the Freeman-Tukey Exact (FTE) inversion. 
Since this formula is applied using an average value of n, 


there may Бе concern that it oscillates between successive 


integer values of n. Figure 3 indicates, using values of 
n=5 and 10, that it is qurte smooth as an interpolaí n 
formula. In the article, Miller suggested using the 


harmonic mean of the n values provided by the time changes. 
This issue is examined later in this chapter. The Miliz 
formula was incorporated in the Marine Corps Inventory 


Attrition Analysis Program provided in Apprendi в. 


De SHRINKAGE VISUALIZATION 

When shrinkage is a fixed amount for all cells, it is 
readily visualized, for all degrees of shrinkage, by the 
linear diagrams presented in the works of Tucker [Ref. 1:p. 
52] and Casella [Ref. 6]. Robinson  [Eef. 2:p. 20] ас 
presented diagrams of the limited translation option and how 
the shrinkage is affected. All of these diagrams refer to 
the transformed scale. 

Of more immediate interest is the question of how to 
visualize shrinkage diagrams of this type in the original 
scale. Illustrations of diagrams of this type are provided 
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using both the basic inversion formula used by both Tucker 
and Robinson, and the FTE formula. Figure 4 shows how the 
diagram by Casella [Ref. бір. 84] looks in the origi romi 
scale using the same value of n for all cells. An 
interesting note to consider is the effect of the difference 
between inversion formulas, 1.е., in the graph πὶ 
combines the Basic and FTE inverses, the FTE inversion tends 
to shrink at a slower rate, as values get further from the 
mean. For a modest transition into the case of non constant 
n, Figure 5 provides for two values of n. Finally, Figure 6 
illustrates the point with many different values of n, taken 
from actual Marine Corps officer inventory data. "Ines 


case, note that Basic inversion shrinks to one original 


scale value while the FTE inverse is not as focused. 


Е. CHOICE OF AVERAGE INVENTORY VALUES OVER TIME 

A Fortran simulation was developed, included as Appendix 
C, to test the performance of each of several common means 
(arithmetic, geometric and harmonic) in the comparison of 
the inverted estimates with the empirical y/n. Values of п 
were produced from a Poisson distribution with rate A (1 sx A 
< 20), and attrition values, y, were produced LL6. 
Binomial (n,p) distribution where a was in the range 01 565 
2.40. The n and y values were then transformed by using the 
Freeman-Tukey formula. An exhaustive study of inversion 


variations with both the FTE formula ama the baste fornu s 
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used by Robinson were compared. Table 3 contains selected 
outputs which are typical of the broader study and shows how 
the means performed against each other in three different 
computations listed below. 


1. The BASIC Method utilizes the Basic inversion formula 
and the average transformed scale value. 


[2 


The ΕΤΕΙ Method utilizes Equation 2.7 with the t 

values equal to the average (over replications, 1.8., 

time) of x/(n + 1/2):/2 and = given by Esquationa wz 

More explicitly, see Equation 2.16. 

3. The FTE2 Method utilize Equation 2.7 with Zen: 

values equal to the average (over replication, i.e.. 

time) of x directly, and x given by качастоп mii 
followed with a division Dy as ze = 

In the context of this study, the arithmetic mean tended 

to perform as well or better than the geometric or harmonic 

means. The choice of arithmetic mean also circumvents the 


problem of an undefined term when ni = 0 15 encountered 


while using the harmonic of geometric means. 


F. QUESTIONS RELATING TO VALIDATION 

The pilot studies by Tucker and Robinson both used four 
years (1977-80) for estimation of the attrition rates and 
the remaining three years (1981-83) for validation. The 
estimation rates were used in an attempt to forecast the 
values for the three validation years. The quality of the 
forecasts deteriorates as lead time increases in the future. 
This was anticipated, in fact, Rowe et al., at NPRDC have 
shown a time series effect in their work forecasting 


attritions in the Navy (Ref. 7] and are currently 


20 


Table 3. 


2.0 


2.0 


2.0 


2.0 


4.0 


4.0 


4.0 


4.0 


8.0 


8.0 


8.0 


8.0 


SIMULATION 


UTILIZATION 
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40 


40 


40 


40 


40 


40 


40 


40 


40 


== —— mm mm i i cm mr mmm a eee _ хш — 


0.10 


0.20 


0.40 


0.05 


0.10 


0.20 


0.40 


0.05 


0.10 


0.20 


0.40 


BASIC 
FTE1 
FIE? 
BASIC 
PTEI 
ETE? 
BASIC 
ETE! 
БЕГЕ? 
ВА51С 
PTEI 
EIE? 
BASIC 
PIEI 
FTE2 
BASIC 
LI LC | 
BLE? 
BASIC 
ΕΤΕΙ 
μα” 
BASIC 
rie 
Ill 
BASIC 
INE 
LI > 
BASIC 
ΕΤΕΙ 
ГОН? 
BASIC 
L ICI 
PTFE? 
BASIC 
FTE1 
ШЕ) 


0.1068 
0.0252 
0.0181 
0.1399 
0.0535 
0.0489 
0.2402 
0.1900 
0.1657 
0.3713 
0.3442 
0.3525 
0.0750 
0.0235 
0.0200 
0.1280 
0.0768 
0.0726 
0.1965 
0.1446 
0.1485 
0.4382 
0.4564 
0.1280 
0.0552 
0.0181 
0.0156 
0.1226 
0.0922 
0.0828 
0.2035 
0.1778 
0.1708 
0.3799 
0.3700 
0.3662 
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0.1176 
0.0159 
0.0205 
0815537 


0.0458, 


0.0558 
0.2629 
0.1829 
0.1869 
0.4038 
0.3403 
Dilo 
0.0815 
0.0200 
0.0221 
0.1387 
0.0724 
0.0798 
022125 
0.1406 
0.1628 
0.4694 
0.4558 
0.3658 
0.0580 
0.0166 
0.0165 
0.1285 
0.0904 
0.0873 
0.2131 
0.1762 
ტი) ს. 
0.3962 
0.3693 
0.3837 


OUTPUT OF MEAN VALUES FOR INVERSION 


0.1308 
0.0122 
0.0245 
0.1708 
0.0369 
0.0648 
0.2906 
0.1739 
0.2144 
0.4429 
058595 
0.4213 


0.0921 


0.0147 
0.0258 
0.1564 
0.0655 
0.0922 
0.2386 
0.1559 
0.1569 
ο ο. 
0.4549 
0,5225 
0.0616 
0.0148 
0.0177 
0.1363 
0.0550 
0.0933 
072295 
0.1741 
0.1915 
0.4170 
0.3684 
0.4065 


recommending a second order auto regression scheme for this 
purpose. In this regard, the work of Tucker and Robinson 
with validation may be viewed as forecasting by persistence, 
and in light of the NPRDC work, can hardly be expected to 
perform well in an absolute sense. Their work does, 
however, provide some comparisons that can be useful in 
choosing among alternative techniques, and much can be 
learned from them. 


The first measure considered is the Mean Squared Error 


(ЛЕР) Ж о forecasts. The shrinkage estimation theory 
provides that these values should be about unity йа 
persistence, independence over time, was agreeable. That 


15, useful validations can be obtained without modeling the 


time dependence of the attrition process. Only for a very 
few cases does this hold e; Warrant Officers, see 
Tucker rn However, if one method of estimation consistently 


produces noticeably smaller MSE values than do the other 
methods, this would provide sufficient support for thas 
method's continued development. 

The second measure is designed to consider performance 


in the original scale andis patterned after a chi-square 


ptg ecc 
k = 2 
= (actuals forecasts) (2.8) 
i Foro l sho 
If the k cells are independent, then ecuation 2.8 15 


approximately a chi-square random variable with k degrees of 


freedom. Hence, the man value is about k and the standard 
deviation is (2k)!/?. Once again, we generally do not see 
values this small, but use it for comparison purposes. As a 
side note, the value of k when using the MLE is generally 
Smaller than that of the other estimators due to the 
omission of those cells with MLE equal to zero over all 
estimation years, i.e., zero forecasts. This ааа апа 
unnatural advantage for the MLE, which must be kept in mind 
when interpreting the results. 

Two disturbing features have emerged from the 
exploratory analysis. Often a technique that performs well 
for some cases using the MSE measure is not comparably 
supported by the chi-square measure. It seems the two 
measures do not weight the common features equally, and 
further study is necessary to develop understanding of the 
results. 

USO with personnel at NPRDC concerning these 
matters has led to the introduction of a third measure, 
which will be called the Mean Absolute Deviation (MAD). 
This value is computed in the original scale. Moreover, 
since the cost structure of over-forecasting is not 
symmetric, we separate the two parts into the average 


overage deviation given by: 


k 
ODEV = 22 [forecasts - actuals]* (2.9) 


1 
κ 


Yon 
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and average underage deviation 


k 
UDEV = x > [actuals - forecasts]}* (2 34290 
1 
where the "+" refers to the positive part, as computed 
separately. Finally, 
MAD = ODEV + UDEV (2.223 


This measure also has shown that in most cases, the 
overage exceeds the underage value, which in turn may lead 
to forecasting the need for an excess number of Marine 


officers. 
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DEE NMPIRIERTFBRIESFESTIMATION 


A. BACKGROUND 

Although the studies of Tucker and Robinson showed 
improvement over the current methods, no clear procedure was 
established and no plans for implementation were drafted. 
Problem areas were identified while working with small cell 
inventories and low loss rates. In the transformed scale, 
the James-Stein assumpticn of normal distribution of Е 
observations with equal variance did not hold well in such 
cases. An alternate method of loss rate estimation is 


required to deal with this special case. 


Bee INTRODUCTION 

Empirical Bayes is simply a generalization of the James- 
Stein estimator for use with small cell inventories and 
unstable variance. Random variables are assumed to be of 


mee form: 


იძია ο i) τπτ 2 E AP 


where ©; are the unknown parameters to be estimated and di 


are known variances. The Bayesian assumption is that the 
©; '5 are themselves from some distribution, in this case 
also from the normal distribution with mean, v, and 
Martance, Pia. The proportionality constant pi 15 


Dermalized so that 


° 
ІІ 

^ 
j 
T 

pa 
ο) 
x 
ω 
ιο 


where k is the number of observations or cells. Thus, 
Q, 7 N(v,pia), τα ο ο... (3. 


This distribution is known as the prior distribution some 
and represents the additional assumed information on the 
value of © before observing X whose distribution depends on 
@ (Ref. 8:p. 553]. The use of this information produces an 
estimate which is a weighted average of the 1. ior mean. v, 


and the sample estimate x. In the empirical Bayes 


Q 


onto oa 
the values of v and a, are also estimated from the data 
since all the information necessary is contained in the 
marginal distribution of Xı (unconditional on 061) [ке k P 
83]. This chapter focuses on that method. 

For a more detailed explanation of the Bayesian family 


of estimators, the reader is referred to Casella [Pef. 6], 


Carter and Rolph [Ref. 4] and Efron and Morris [Ref. 2]. 


C. THE ESTIMATION METHOD 
Fay and Herriot [Ref. 10] discuss the empirical Bayes 
technique and provide an estimate for the prior mean as 
> 2 / 3“. c 


у^ = 1 нае (3.49. 
E 1/(a* di) 


1 
where a* is an unbiased estimate of a. The empirical Bayes 
estimate of ©, commonly called the mean of the posterior 


26 


where a* is an unbiased estimate of a. The empirical Baves 
estimate of O, commonly called the mean of the posterior 


distribution becomes 
a ას. # 
πον. var «απ ον. 


Mims 16 the mean of the conditional distribution of 9 given 
еее. 8:p. 556]. The weights of each of the prior and 
sample means are determined by the prior and observed 
variancss. 

As a measure of the worth of this estimator, when the 
number of x's is 2 4 [Ref. 11], the risk, or sum of squared 
differences between the true unknown parameter and the 
empirical Bayes estimate, is less than that using the 
observed x's, for all 01. Efron and Morris [Ref. 9] provide 
r iorous proof. 

An important result of the empirical Bayes estimator is 
that shrinkage values are no longer constant, but depend 
ШЕСІ intormation contained in each cell. Additionally, the 
point toward which the cell values are shrunk is the 
weighted mean, not the grand mean of the observations. 
Appendix A provides the derivation of the empirical Bayes 
estimator, variable shrinkage factors and calculation of 
the prior variance used in this study as they pertain to the 


Marine Corps officer attrition rates. 
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IV. VALIDATION ომლის. ТЕЗОЮ 


A. GENERAL 

The outputs displayed in this chapter are excerpts from 
that produced by the Marine Corps Personnel Inventory 
Attrition Analysis (MCIAAP) (see Appendix B). Initially, 
comparisons will be made with the output displayed by 
Robinson to identify the validity of the Fortran code with 
the APL code used by Robinson, then new approaches will be 
displayed in an attempt to improve the results. The Fortran 
code in the MCIAAP incorporates most of the procedures 
studied by Tucker and Robinson. Specifically, the MLE, 
Transformed Scale Cell Average (TSCA) and the James-Stein 
(JS) estimators, shrinkage procedures and a basic inversion 
formula. The FTE inversion formula, the new measures of 
effectiveness from NFRDC, are also included. For the MAD 
measure of effectiveness, only maximum likelihood and the 
modified empirical  Bayes variables will be presented, since 
these showed the most promise over all. All other variables 
are displayed by the MCIAAP if interest is warranted by the 
reader. Chapter II describes all of these methods in some 


detail. 


Be: HISTORICAL COMPARISONS 
Tables 4-9 were developed using the same aggregation 
scheme and format as Robinson [Ref. 2:pp. 39-41), altheien 
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not all of the variables are included (MLE, TSCA. and James 
meein are the only variables duplicated in this paper). 
Additionally, Robinson presented the number of cells, ask 
values, used for evaluation of each estimate in the 
aggregate Figure of Merit (FOM) display. These k values are 
not repeated in the current FOM tables, but are included, as 
parenthetical values, and discussed in the new MAD measure 
of effectiveness tables. Since this is the case here, the 
reader is reminded that the apparent competitive performance 
of the MLE in the following cases is due, in part, to the 
possible advantage of having a smaller number of cells used 
TRIS MOE's computation. 


Since the computations were identical for both studies 


for these tables, the output was expected to be nearly 
identical or some pattern discernable between the two. In 
very few cases was this found to be true, and after 


verifying that the algorithms coded both in the APL and 
Fortran versions Were correct, we have no reasonable 
explanation for the discrepancies of the programs. It is 
suspected that the transfer of raw data into an APL 
workspace, accomplished by Tucker and Robinson, may be at 
EEJance with that into Fortran files, done by contractor. 
1. Aviation Aggregate Case 

In the case of aviation FOM's, displayed in Table 4, 

1981 values for LtCol's in the transformed scale are 


πο. ο σαι to Robinson's [Ref. 2:p. 39) while the last two 
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Table 4. 


AVIATION FIGURES OF MERIT 


Ist Lt 


MLE 
TSCA 

JS 

MOD JS 
EB 

MOD EB 


Lt Col 


MLE 
TSCA 


MOD EB 


Ist Lt 


MLE 
TSCA 
TSCAM 
JS 

JSM 

MOD JS 
MOD JSM 
EB 

EBM 
MOD EB 
MOD EBM 


Lt Col 


МЕР 
ISCA 
TSCAM 
JS 

JSM 

MOD JS 
MOD JSM 
EB 

EBM 
MOD EB 
MOD EBM 


1981 


[082 


1983 


TRANSFORMED FOM 


4.0563 
3.7420 
3.9694 
4,7352 
4.0753 
3.1533 


.3664 
.7168 
.7 369 
.7208 
8482 
‚3501 


L, Q tn £. 


, 


21.0501 
56.4603 
53.7884 
58.0150 
55.3844 
63.0345 
60.6413 
28.8824 
ΠΠ 
19.7653 
20.7987 


55 73515 
38.4478 
35.0390 
37.8248 
33.2016 
38.0291 
34.5466 
40.3421 
36.7498 
40.2056 
185.7476 
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10.6590 
10.3636 
10.5671 
11.2615 
10.7015 

5.0568 


8.5033 
10.2994 
10.2468 
10.1865 
10.2293 

3.9824 


ORIGINAL FOM 


43.1546 
95.7510 
90.5164 
95.2150 
91.9406 
99.5419 
96.5615 
96.2449 
9219699 
38.4637 
36.3673 


59.4200 
57.0937 
48.9374 
3251 15 
46.9901 
9173225 
45.1554 
54.9389 
46.9657 
70.5331 
65,1553 


11.1338 
10.5711 
11.1626 
12.1206 
607 

6.2870 


8.1506 
1023355 
10.2972 
10.2756 
10.3554 

2.7186 


53.0528 
98.4715 
955920 
100.1664 
97.0163 
105.4897 
102.5172 
101.2482 
98.0958 
49.7393 
48.1659 


22.1156 
69.0641 
ο 
67.1696 
027129 
63.3908 
οσο. 
67.7091 
3789307 
60.5749 
168.3976 


years are smaller. 


years. The modified empirical 


best overall yet still may be 


E: bé close to unity). TEI IS 


most cases, the risk values 


forecasts, 
Шггвсасе with time. 


πατε.» original scale, 


Robinson's for all years, 


were all larger. The 


variables except for 


validation years. MLE 


and are approximately within 


deviation of the number 


LtCol's, FTE inverted JS was 


БЕ ТЕРЕ Ρο Ss are larger for all 


yet the Mod EB FOM decreases as 


MLE tended to be smaller 
while TSCA 
ΕΤΕ inversion 
Mod EB 


and Mod 


of cells 


best 


chree 


Bayes (Mod EB) values were 


too large (ideal values tend 
"ესო. ა-ა ილ со Note that in 
tend to increase with future 


future forecasts 


than 
and James-Stein (JS) 


и 


worked best 


a. 


Tho thue 1 IL an le En 


EB worked best for 1stLr's 


the range of one standard 


Im tlhe computation. Fer 


with MLE a competitive 


second, yet all were not in the desired range. 


In the MAD measure, 


between future forecasts 


the two ranks (FTE Mod EB was 


Pest LtCol's). The mix in 


BEccorable in the first year for all cases, 


Sells (К) 


Sell numbers in most cases, 


the performance of these 


Suspect. The second year 
displays a preferred ratio 
value may still be too large. 


რე. 


Table 


(no set pattern) 


overage 


estimates 
values 


of K 


5, results are mixed 


and variables £or 
best for 1stLt's and MLE was 
is most 


and underage 


yet the number of 


forecasted as overage exceed the underage forecast 


Dis iS eause for concern and 


іп this aggregate is 


for LtCol estimates 


Values, Dut the total MOE 


Table 5. MEAN ABSOLUTE DEVIATION (AVIATION) 


1981 1982 1983 
Ist rt 
MLE 
UNDERAGE 0.1256( 3) 0.0439( 1) 0.3069( 1) 
OVERAGE 2.3985( 8) 3.9267(10) — 4.2522(10) 
TOTAL 2.5240(11) — 3.9706(1) 4.359111), 
MOD ED Р | 
UNDERAGE 0.1822( 4) 0.0410( 1) 0.3281( 1) 
OVERAGE 2.3625( 8) 3.9040(11) 4.2351(11) 
TOTAL 2543702) 3.9451(12)  4.5633(12) 
MOD EBM 
UNDERAGE 0.2518( 5) 0.0580( 2) 0.3420( 1) 
OVERAGE 2.1654 7) 3.6427(10) 3.9677(11) 
TOTAL 2.417212) — 3.700712) 4.3098(12) 
Lt Col 
MLE 
UNDERAGE 0.6005( 5) 1.5139(10) 0.0396( 1) 
ОУЕВАСЕ 0.7925( 8) 0.3928( 3) 1.3467(12) 
TOTAL 1.3930(13) — 1.90673) — 1.386313) 
MOD EB 
UNDERAGE 0.6053( 6) 1.4717 9) 0.0841( 2) 
OVERAGE 1.3083(10)  1.0970( 8) 2.2331(14) 
TOTAL 1.9136(16) 2.5687(17) 32317316) 
МОР ЕВМ 
UNDERAGE 0.6890( 6) 1.5366(10) 0.1167) 2( 
OVERAGE 1.0120(10) 0.7179( 7) 1.750114) 
TOTAL 1.7010(15) 1.8669( 16) 


2.2546(17) 


2. Combat Support Aggregate Case 


Table 6 displays the combat support aggregate output 
from the program. In the transformed scale comparison with 
Robinson [Ref. 2:p. 40], the values for the two ranks take 


opposite directions. Where the 1stLt FOM values are all 
‚Ф 
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Table 6. 


COMBAT SUPPORT FIGURES OF MERIT 


MOD EB 


Lt Col 


MLE 
TSCA 


MOD EB 


st Lt 


MLE 
TSCA 
TSCAM 
JS 

JSM 

MOD JS 
MOD JSM 
EB 

EBM 
MOD EB 
MOD EBM 


Lt Col 


MLE 
TSCA 
TSCAM 
JS 

JSM 
MODJS 
MOD JSM 
EB 

EBM 
MOD EB 
MOD EBM 


1981 


1982 


1983 


TRANSFORMED FOM 


21977 
1.6329 
1.6925 
2.3684 
2:3323 
2.2293 


1.3103 
0.7502 
0.6421 
0.7009 
0.7736 
1.5691 


111.4260 
73.5008 
79,0947 
76.9945 

429.8503 

143.1700 

168.9585 

10229525 

108.3578 
65.7251 

176.3961 


36.0152 
61.7769 
059525 
61.0511 
66.3509 
66.5470 
51.4722 
69.7508 
Ἡ] 5458 
62.6316 
83.1337 
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ORIGINAL FOM 


112.1917 
90.0523 
95.8244 
89.8064 
84.4175 

113.0083 

114.7083 

119.0900 

l11.2701 
63.9054 
395700 


66.0261 
52.0346 
49.6307 
49.2280 
47.8490 
49.9366 
30.2657 
59.9812 
50.0281 
501.5 
4S.6576 


70.2767 
73.9380 
80.9764 
"0 
65.2639 
90.0783 
920399 
ο 
36.5133 
45.2719 
43.0037 


42.4097 
TOL, 
107.1016 
96.7861 
118.2564 
99.2566 
105.5460 
104.5311 
109.2012 
96.8693 


127.1541 


larger (except for MLE ain 1991) Пе LtCcol νετ’. 28 


(i 


ES 


smaller. Overall, ail estimates are competitive with eac 
other over all validation years with TSCA performing best 
for the 1stLt's and James-Stein outperforming the rest in 
the LtCol's. Note here also that the Mod EB estimate did 
not fair well in either case, and was the worst estimate in 
risk foriti  πσο с. 

In the original scale , a different and more random 


effect is noticed between Table 6 and the results of 


= 


pu 
- 


Robinson. Risk values tended to fluctuate with for 


iD 


U) 


lead time between estimates for 1stLt's where the LCCC1 
values were more stable, yet were quite different from the 
transformed scale pattern. In this case, Mod EB performed 
well for  1stLt's, and was the only one within the desired 
range of standard deviations. MLE for the LtCol's was also 
the only estimate desirable fcr this measure. The Miller 
inverse showed poorly in this aggregate for all estimation 
methods, in every case exceeding it's counterpart. 

In the MAD measure of effectiveness, Table 7, MLE 
performed the best in both ranks. The Mod EB effect is 
competitive in +he 1stLt's case, and all variables in the 
Ltcol case are closely grouped. The number of cells are 
heavily weighted towards an overage forecast in all three 
validation years, throughout this aggregate. Once again, 


this is not a desirable quality for a MOE. 
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Table 7. MEAN ABSOLUTE DEVIATION (COMBAT SUPPORT) 


Ist Lt 
MLE 


UNDERAGE 
OVERAGE 
TOTAL 


MOD EB 


UNDERAGE 
OVERAGE 
TOTAL 

MOD EBM 
UNDERAGE 
OVERAGE 
TOTAL 

Lt Col 

MLE 
UNDERAGE 
OVERAGE 
TOTAL 

MOD EB 


ENDERAGE 
OVERAGE 
TOTAL 


MOD EBM 


UNDERAGE 
OVERAGE 
TOTAL 


0.3737(14) 
0.6731(19) 
1.0468(33) 


0.4028(15) 
0.7311(40) 
1.1339(55) 


0.4592(16) 
0.5930(39) 
1.0522(55) 


360( 9) 
510(14) 
$70(23) 


0.1 
0.1 
02 


0.1638(11) 
0.2232(29) 
0.3870(40) 
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0.1617( 8) 
0.9574(25) 
1.1191(33) 


0.1275( 9) 
0.9972(45) 
1.1246(54) 


0.1603( 9) 
0.5318(45) 
0.9921(54) 


0.1688( 7) 
0.2591(16) 
0.4280(23) 


0.1901( 7) 
0.2240(35) 
0.4142(42) 


0.1781( 8) 
0.8835(25) 
1.0616(33) 


0.1427( 9) 
0.9422(38) 
1.0849(47) 


0.1856(10) 
0.7575(37) 
0.943 1(47) 


0.05-1( 5) 
0.3415(18) 
0.3957(23) 


0.1162( 6) 
0.4548(29) 
0.5710(35) 


0.1408( 6) 
0.2589(29) 
0.3996(35) 





CJ 


Ground Combat Aggregat=] ase 


In the final aggregate studied by Robinson [Ref. 


2:p. 41], that of ground combat MOS Ss, the results again Е 
displayed in Table 8, were mixed between computational 
methods. In the original scale, value differences tended tc 


be larger than Robinson's in all cases for lstLt's sieve 
a preporticnal growth over Επ Conversely, LtCol values 
were smaller іп the first two validation years and slightly 


larger in the final year. τους Pooks good for 


rh 


with the others being competitive, yet the size of the 
values may nct be desirable. Mod EB turned out the best for 
LtCol's, but all other estimates are competitive and in the 
desired range. 

On the original scale, all values for the estimates 
for lstLt's are too large, with Mod EB performing Welle 
that regard. LtCol values are all larger than Pobinson's, 
but the estimates are all within range of one standard 
deviation and, therefore, competitive. MLE, once again, 
shows well for this rank. In opposing fashion to that m 
the combat support aggregate, the FTE inversion outperformed 
all Basic inversions for all variables. 

The MAD measure of effectiveness, Table 9 confirms 
the results found above ın that Mod EB was best for 1stLt's 
and. MLE азай е1 γος ΠΟΘ ss The values for 1stLt were 
found to be too large overall, while all estimates for ЕСО 


were competitive and in the desired range. The mix between 
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MOD EB 


Lt Col 


MLE 
TSCA 


1st Lt 


MLE 
TSCA 
TSCAM 
JS 

JSM 
MOD JS 
MOD JSM 
EB 

EBM 
MOD EB 
MOD EBM 


Lt Col 


MLE 
TSCA 
TSCAM 
JS 

JSM 
MODJS 
MOD JSM 
EB 

EBM 
MOD EB 
MOD EBM 


1981 


Table 8. GROUND COMBAT FIGURES OF MERIT 


1552 


TRANSFORMED FOM 


2.5764 
1.9741 
21737 
3.0859 
2.5558 
3.1538 


87.3362 
S7 I912 
83.7648 
95.4717 
92.0238 
124.6490 
134.9472 
109.6349 
105.5708 
64.5647 
63.6623 


35.0538 
55.9649 
43.5234 
56.3346 
43.9957 
60.0730 
52.5940 
62.7441 
50.4049 
48.0871 
29 1231 
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4.1136 
3.8072 
3.9738 
4.7919 
4.3026 
3.3202 


1.9674 
ο 
22119 
2351 
2.2464 
2.3930 


ORIGINAL FOM 


119.2637 
116.6767 
ΠΠ» 
122.9425 
117.5264 
147.6611 
142.5461 
150952] 
131.0400 

81.2071 

77.8159 


45.8172 
198.8659 
201.9953 
198.9037 
2013299 
202 2959 
210 7551 
205.4249 
209.3533 
197.7316 
323.5654 


6.8744 
ΠΠ; 
6.3970 
7.5619 
6.9607 
4.7069 


G2 us‏ ړن دی دی رح 
GU LGI) — I. Un —‏ 
[S oo ο ο οὐ‏ 
O Un UY Un DE‏ 


CO თ) C5 Us Ln O 


143.9128 
145.6236 
142.0902 
(52554 
149.3965 
150-3507 
178.6173 
102528] 
165.S1S9 
101.9511 
100.7150 


53.0111 
123.1306 
122.8648 
1222750 
115.7143 
129-9299 
121.4254 
126.6544 
125.7642 
Π1 555 
242.4398 


over-estimates and 


cells  forecasted 


desirable one, 


validation years. 


preceding aggregates. 


or 


the same 


This also 


under 


under-estimates, ige, 


the 


the actuals, is 


mix holds 


was the 


true 


finding 


Table 9. MEAN ABSOLUTE DEVIATION (GROUND COMBAT) 


lst Tt 
MLE 


UNDERAGE 
OVERAGE 
TOTAL 


MOD EB 


UNDERAGE 
OVERAGE 
TOTAL 


MOD EBM 


UNDERAGE 
OVERAGE 
TOTAL 


Lt Col 

MLE 
UNDERAGE 
ОУЕКАСЕ 
TOTAL 

MOD EB 
UNDERAGE 
OVERAGE 
TOTAL 

MOD EBM 
UNDERAGE 
OVERAGE 
TOTAL 


1981 


0.1525( 5) 
1.9322(30) 
2.0847(35) 


0.2340(10) 
1.7743(33) 
2.0083(43) 


0.2618(10) 
1.6492(33) 
1.9111(43) 


0.2313(10) 
0.4182(25) 
0.6495(35) 


0.2315(11) 
0.6392(34) 
0.8706(45) 


0.2650(12) 
0.4071(33) 
0.6720(45) 
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1982 


0.1855(10) 
2.5979(25) 
2.7834(35) 


0.1323( 8) 
2.3086(32) 
2.4409(40) 


0.1705(10) 
2.1738(30) 
2.3443(40) 


0.3034(11) 
0.4202(24) 
0.7236(35) 


0.3790(15) 
0.4427(31) 
0.8217(46) 


1983 


0.2466(10) 
3.1691(25) 
3.1157(35) 


0.1368(10) 
0.6453(25) 
0.7821(35) 


0.1426(12) 
0.9525(31) 
1.0950(43) 


0.1864(12) 
0.6535(31) 
0.8399(43) 


number of 


tor 


I. Test Aggregate Case 
In an effort to experiment with aggregate groupings, 
another case was built upon which all estimation methods 
were applied, and the results compared with previously 
studied groups. The grade was held fixed, in this case the 
ο οσο Captain was used, with a LOS in the eighth year. 
Several MOS's were chosen so that to insure an aggregate 


Wien non-zero inventories. In the transformed scale, Table 


10, all estimates have small and competitive values, with 


М 


ЕВ and MLE doing well. In the original scale, onc: 


( 


agaln, all estimates are well within range, but Mod EB and 


MLE were the standouts. ე: . ілуегесоп ортоп ата пос 
fair well with this aggregate. All estimates in the MAD 
measure of effectiveness, Table 11, are very competitive. 


and the first validation year mix (overage and underage) may 
be the best yet. It is interesting to note that the cell 
split in the first validation year between overage and 
ШЧвгапе forecasts is contrary to that of the historical 
displays above. Although the total number of cells used 


(k=6) 15 much smaller than the previous table, it is felt 


(D 


that the aggregate scheme was the major factor in this 


finding, and the results are promising. 


Sym 


Table 10. TEST AGGREGATE FIGURES OF MERIT 


1981 1982 1983 
TRANSFORMED FOM 
Captain | 
МЕЕ 0.9603 2.3463 2.0842 
TSCA 1.4151 3.3033 29152 
JS 1.3753 3.3969 25205 
MOD JS 1.3943 3.S165 2.8829 
EB 1.6840 4.2698 3.2957 
MOD EB 0.8249 |» 1.2950 
ORIGINAL FOM 
Captain s 

MLE 3.0402 4.8104 4.9399 
TSCA 4.6161 10.2779 7.5714 
TSCAM 8.1467 7.6241 8.8337 
JS 4.3970 10.3000 73726 
JSM 7.3199 7.6666 5.3560 
MOD JS 4.2990 10.76-0 7.6371 
MOD JSM 6.0079 5.5070 7.9259 
ЕВ 5.3660 12.0537 8.8679 
EBM 6.7181 9.3915 8.7959 
MOD EB 1.4463 6.0323 3.9066 
MOD EBM 4.6465 4.3554 59173 


Table 11. MEAN ABSOLUTE DEVIATION (TEST CASE) 


Captain 
NILE 
UNDERAGE 
OVERAGE 
TOTAL 


MOD EB 
UNDERAGE 
OVERAGE 
TOTAL 

MOD EBM 


UNDERAGE 
OVERAGE 
TOTAL 


1981 


0.1324) 
0.2575( 2) 
0.3918( 6) 


0.3712( 4) 
0.0852( 2) 
0.4563( 6) 
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1982 


0.0463( 1) 
0.7374( 5) 
0.7837 6) 


0.0000( 0) 
0.9601 6) 
0.9601( 6) 


0.0809 1) 
0.6148( 5) 
0.6957( 6) 


1955 


0.2097( 2) 
0.4391( 4) 
0.6488( 6) 


0.0965( 1) 
0.6306( 5) 
0.7270( 6) 


У. CONCLUSIONS AND RECOMMENDATIONS 


AS > US 

Omlteing Ene fact of differences in computational 
methods between this paper and the Robinson thesis, it is 
clear that the results are still inconclusive at this time. 
In the historical comparisons, the estimates studied varied 
in their performance from one aggregate to another. The 
different risk calculations, in the transformed and origi emi 
Scales, as well as the different measures of effectiveness 
showed little agreement in estimation methods. Long range 
estimation still remains a problem since with time, the 
estimate risks increase significantly. Much aS s sN was 
put on the modified empirical Bayes estimator in this paper 
and those specific results look promising under росе 
different and more stable aggregate conditions. 

One area of mention here concerning the results has to 


do with the composition of inventories used in the selected 


aggregates. It was noticed that the range of inventory 
(smallest to largest value in any particular MOS) was quite 
large, and may be affecting the results. Im fact, Carter 


and Rolph [Ref. 4:p. 3883] managed this point by further 
disaggregating into levels of activity as it concerned their 


work with fire alarm -<tc πα тол? 
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Be CONCLUSIONS 
The results of the previous feasibility studies and this 
paper leads to several conclusions. The chosen aggregate 


still seems to be the key element in the work of attrition 


estimation. It does seem clear that the aggregation scheme 
Should be focused on fixed ranks, in narrow LOS ranges and 
multiple MOS's. This conclusion is based upon the relative 


Si@lecess of the test aggregate utilized in this work, and 
upon Similar success by Carter and Rolph. Additionally, the 
Meet rsioning of inventory levels in the aggregate should 
stabilize the estimation method performance. All the 
estimation method studied to date have their strengths, but 


still, no one method of estimation excels over the others. 


Es RECOMMENDATIONS 
Mo estimation method studied in this paper is 
recommended for implementation at this time. Tiada) τομ, 
no method is recommended for removal from consideration 
since basic aggregation problems are felt to be the largest 
contributory factor to the mixed results found here. These 
recommendations are presented for further study. 
n. The aggregation method offered by Captain Larsen 


[Ref. 3) should be examined further with its 
relevance to the scheme in the test case used here. 


L2 


Utilization of the detailed data provided by NPRDC as 
soon as possible in order to expand the available 
information and open new avenues of aggregaticn. 


e. Re-examine the inventory mean value used in the FTE 
inverse CC determine its future role in the 
estimation process. 
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APPENDIX A. ΕΘΙΜΟ σος 


A. NOTATION 

The following notation is a generalization cf that used 
by Robinson [Ref. 2:p. 71] as it pertains to the aggregation 
scheme used in this paper. The indexing system has been 
changed to reflect the use of officer grades to identify the 


cells. 


|5- 
il 


index of MOS cells in the aggregate. 

J = index of LOS celis in the aggregate. 

k = index of Grade cells in the aggregate. 

t = index of time periods in the aggregate. 
invli,j,k,€) = inventory with mos 1. LOS 3 and Grace 


during year t. 


y(i,j,k,t) = number of attritions in cell (i,j F) dur m 
year t. 
n(i,j,k t) = maximum 209995 "'""” 


inv (E E 
D = matrix identifying cells with no inventory over all 


estimatior years. 


თს) к) O, 1£ cell -is a stru о. 


თე. 1, ТЕ cell is not a Structural zero. 
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B. PREEMAN TUKEY EXACT INVERSION ALGORITHM 
wnis algorithm provides a formula for the inverse of the 
Freeman-Tukey double arcsine transformation to express the 
means of the transformed values into proportions on the 
original scale. 
STEP 1. Transform the data using the Freeman-Tukey 


_ Те агссіпе transformation. 


МТ kt) = დაუ... ი... DERI -1] 


= > 7 
DER Lo 1. ა.ს... -111 ( 


25 = _ 1 - 4 \ 
ο οσα ЕЕРЕЕ ум 
STEP 2. Calculate the average cell inventory and 
transformed values over all Т. 
x(i,j,k) = (1/T)E x(1,3,k), (A.2) 
ДЕ 
ποσο απ ο τι (II) > πη. Πλ (A.3) 
E 
STER Check for inversion values less than zero 
or greater than one. The algorithm is terminated at this 
БЕШПС lf these conditions are met. 
tm = x(i,j,k)/(inv(i,j,k)+0.5)1/2 + п/2, (A.4) 
emp = sim: (1/(1ny(1,3,k) + 1))t/2, (А.5) 
IC. tm < LX), FTE = 0.0, 
IF tm » (n-tmi), ETE S 1.0, 
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STEP 4. Compute the 5-0 ელ... 


If Cos (emi СО 59 - ο 


If Cos წი) Qo зап = -1.0. 


STEP 5. Computer the охе сеч тае сите 
transformed variable 


FTE = 0.5(1-soalas [sin ita) > Әрекет NE 
по. к. (р 


er TRANSFORM SCALE VARIANCE ALGORITHM 

This algorithm produces the variance of the transformed 
data within the unstable variance range where the assumption 
of normality does not hold. As discussed in Chapter II, 
this range o£ instability was determined to be between 1.001 
and 2.2. The method for computing the variance within 0 
range is presented here. 

STEP 1. Compute the transformed scale values to be 
used for the variance. Since x(i,j,k) was computed in the 
FTE inversion algorithm of this appendix, the formulae will 
not be repeated. 


x(t) = > О. Ем (ბ... 
13k 


a 
> 


аа 


4 





= II. ӨЛЕ = LII/ >) · (A.8) 


к 
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STEP 2. Compute the variance of 


values. The formula is derived 


regression equation. The range of the variance 


νι. 05 < var(i,3j,k) s 1.0. 
ЕНІП) = a2z(igg,k)P8(z (i,j,k) = 1)9?, 
MOIST E COD 
where a = 1.6835, Ь1 = -0.8934, апа b2 = 0.8991. 


OI IPTPETICAL BAYES SHRINKAGE ALGORITHM 


This algorithm computes the empirical Baves 


values and estimator of attrition rates. Prior 


this algorithm, 


დას Lhitialization. 
аі - а. 
STEP 2. Define the a and 7 parameters 


empirical Bayes context. 


u = 2 l/(a + var(i,j,k)), 


ijk 
О 


where varí(i,j,k) is defined by formula 


э О ЖМ (ОО A E). 
11 ijk 


y = 


47 


I) თ 


the value for a has been set to 0.0. 


transrormed celi 
Ша Епарвер ЕТ 45 arlınear 


із limited 


(p 


to entering 


(A.109) 


in the 


г> 
|-> 
pa 


LO 


STEP 3. Compute the empirical Bayes prior variance. 


К > ое: 


a 11L | | (A.13) 
ο ο ο ο) 


where ze Y XXII, А MENT 
ijk iik 


STEP 4. Check for the following conditions ss 


branch accord Phnom 


If, a < 0.0 set a = 0.0 and ας со сс -. 
τα, |а H а > “000. return to step 1. 
STEP 5. Compute the empirical Bayes estimator. 
xt?(i,j,k) s e jJx(i.j.X) + (κ... 


atvar (1, 1. K) 


[ var (i jJ K) a 


atvar(i,j,k) 


STEP 6. Compute the empirical Bayes shrinkage set. 


x` 


shrtB(1,j,k) = Var gd: NUMEN cM (ბ... 


E MODIFIED EMPIRICAL BAYES SHRINKAGE ALGORITHM 
This algorithm differs only slightly from that oom 


empirical Bayes shrinkage algorithm listed above, and for 
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that reason, only the pertinent steps will be displayed 
here. 

Consider a particular cell and its transformed value 
cor t=1,...,T in the estimation set. Me, form a 


weighted time average, set 


απ ο ο ο inv(t)]t/? е 


This leads to a modification of step 1 in the following 


Form. 
STEP 1. Compute the weighted transformed scale 
average. 
ces 1 
xt =n Do MEN ს... 
ү 


Since the weights will modify the variance of the weighted 


transformed scale values. 


თა ა = Var( 1, j,k) /(02.5 + inv(t)) (А.18) 
var [xt (t) ] = SE Asya tE = УБИТ (A.19) 
t 


The extreme right of formula A.19 serves to define vt, 
eye approximation for the weighted varlance in the 
meaeadation set for that cell. 

The values of xt and var(xt) are carried cehroughout tnis 
version of empirical Bayes estimation to form the values 


XTEB for the empirical Bayes attrition rate. For the 


purposes of inversion to the original scale, the usual Basic 
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formula is modified since the weschting factor, 0 S MIU 
is already incorporated into the transform value. The new 
Basic inversion for this version is simply as follows: 

р = 1/2[1 + sin(XTEB)] (А.200 
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APPENDIX B. MARINE CORPS PERSONNEL INVENTORY ` 
ATTRITION ANALYSIS PROGRANI 


PROGRAM MARCOR MAR00010 
jeyeyesedede desee je dee dfe je de ενος ΛΑΟ ΑΚ ΚΑΛΑ ΑΚΟΑ ΚΑΚΟ ΚΟΚΚΟΙ ΚΚ ΧΑΑ  ΜΑΒΟΟΟΖ2Ο 


* MARINE CORPS PERSONNEL INVENTORY ANALYSIS PROGRAM * MARO0030 
* *  MAROOO040 
* PURPOSE: TO PROVIDE FUTURE PERSONNEL LOSS ESTIMATION OF THE * МАК00050 
* MARINE CORPS OFFICER STRUCTURE BY SEVERAL DIFFERENT  *  MAROO060 
* ESTIMATION SCHEMES. | * MAR00070 
E PROGRAMMERS: LUIS URIBE, INDEPENDENT CONTRACTOR * МАК00080 
* CAPT C R DICKINSON, USMC * МАК00090 
* DESCRIPTION: WRITTEN IN FORTRAN 77 FOR THE ІВМ 3033 MAINFRAME * MAROO100 
* COMPUTER RESIDENT AT NPS, THIS PROGRAM INCORPORATES % МАК00110 
* MET THODSTEXELORED BY MAJOR D.D. TUCKER , MAJOR J.R. * МАК00120 
* ROBINSON AND CAPTAIN DICKINSON. ESTIMATION SCHEMES * MAR00130 
* INCLUDED ARE: * MAROO140 
* MAXIMUM LIKELIHOOD ESTIMATION (MLE) * МАК00150 
* TRANSFORMED SCALE CELL AVERAGE (TSCÀ) * МАКОО160 
* JAMES -STEIN *  MAROO170 
* LIMITED TRANSLATION JAMES-STEIN *  MAROO180 
* EMPIRICAL BAYES ESTIMATION * МАК00190 
* RISK ANALYSIS WAS PERFORMED IN BOTH THE ORIGINAL AND *  MAROO200 
* TRANSFORMED SCALES. ON THE TRANSFORMED SCALE, * МАК00210 
* ACHIEVED BY USING THE FREEMAN-TUKEY DOUBLE ARCSINE % MAROO22O 
* TRANSFORMATION, THE RISK WAS DEFINED BY THE AVERAGE  *  MAROO230 
* SQUARED DEVIATION OF ACTUAL AND ESTIMATED VALUES. * МАКОО240 
* FOR THE ORIGINAL SCALE, THE TRANSFORMED VALUES WERE  * МАК00250 
* FIRST INVERTED USING TWO TECHNIQUES. AN AD-HOC * MAROO260 
* FORMULA USED BY TUCKER AND ROBINSON AND AN EXACT * МАК00270 
* INVERSE FORMULA FROM JOHN J MILLER. TWO MEASURES OF * МАК00280 
= EFFECTIVENESS, CHI-SQUARE GOODNESS OF FIT AND MEAN თ МАКОО290 
* ABSOLUTE DEVIATION (MAD), WERE THEN EVALUATED. * МАК00300 
* RESULTS ARE PRINTED, WITH ALL ARRAYS COMPUTED, AT THE * МАК00310 
* COMPLETION OF THE PROGRAM. * МАКООЗ20 
* VARIABLE DESCRIPTION: * МАКООЗЗО 
* THE FOLLOWING LIST OF VARIABLES, WITH A BRIEF * MAR00340 
* DESCRIPTION, FORM THE BASIS FOR THE VARIABLE NAMES * МАК00350 
* FOR THE ENTIRE PROGRAM. * MAROO360 
* MLE - MAXIMUM LIKELIHOOD ESTIMATOR * МАКООЗЈО 
* TSCA - TRANSFORMED SCALE CELL AVERAGE * МАК00380 
* J - JAMES-STEIN ESTIMATOR *  MAROO390 
* JILCEUAMESSSDEINOESTIMATOR MODIFIED BY * МАК00400 
© A VARIATION ON SHRINKAGE % МАК00410 
* JLT - LIMITED TRANSLATION JAMES-STEIN * МАК00420 
* ЕВ - EMPIRICAL BAYES * МАК00430 
* TEB - EMPIRICAL BAYES MODIFIED BY A * MARO0440 
E VARIATION ON % МАК00450 
ж TRANSFORMED SPACE VARIABLES * MAROO460 
ж THE FOLLOWING ARE UNIQUE VARIABLES IN THE TRANSFORMED * MARO0470 
* SPACE: * MARO0480 
E X - VALUE OF THE FREEMAN-TUKEY * МАК00490 


8 
oT 
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ж cb ob ob ob cb cb cb cb cb GB Gb ob ob ob cb cb cb ob ob OR Gb Gb cb cb ob OR cb cb cb cb cb Ob cR cb ob OR Gb OR Gb Gb Gb GB GB GR Gb HH AH HH πα 


TRANSFORMATION 


XB > 


(Х-ВАК 


) AVERAGE TRANSFORM VALUE 


OVER ESTIMATION YEARS 
XBB - GRAND MEAN OF TRANSFORMED VALUES 


SSE 2 
SSI =- 


SSB - 


SUM OF 
(X AND 
SUM OF 
(X AND 
SUM OF 


SHR - 


SHRI 


SHRED 
SHRTEB - 


SQUARED DIFFERENCES, ERROR 
X-BAR) 

SQUARED DIFFERENCES, TOTAL 

XBB) 

SQUARED DIFFERENCES, BETWEEN 
(SST-SSE) 

JAMES-STEIN SHRINKAGE VALUES 
MODIFIED JAMES-STEIN SHRINKAGE 
VALUES 

EMPIRICAL BAYES SHRINKAGE VALUES 
MODIFIED EMPIRICAL BAYES SHRINKAGE 
VALUES 


THE FOLLOWING ARE MODIFIED BASIC VARIABLES IN THE 


TRANSFORMED SPACE. THIS IS DONE BY PLACING AN 'X' FOR 


A PURE TRANSFORMED VARIABLE, OR AN 'R' FOR A RISK 
VARIABLE IN FRONT OF THE BASIC VARIABLE. 


MLE 
TSCA 
XJ 
И! 
ХЈІТ 
ХЕВ 
ХТЕВ 
ORIGINAL SPACE VAR 


IABLES 


RML (MLE) 
RTSCA 

RJ 

RODI 

RJLT 
RXEB 

RTEB 


THE FOLLOWING ARE MODIFIED BASIC VARIABLES IN THE 
ORIGINAL SPACE. SIMILAR TO WHAT IS DONE IN THE 
TRANSFORMED SPACE, A 'P' SIGNIFIES A PURE INVERTED 


VARIABLE, AND AN 
'M' AT THE END OF A VARIABLE NAME IDENTIFIES A MILLER 


I 


INVERSE VARIABLE. 


PMLE 
PJ 
PJH 
p 
PJ1M 
PJLT 
PEB 
PEBM 
PTEB 


PTEBM 
THE FOLLOWING LIST IS OF PREFIXES USED TO DESCRIBE 
THE MEAN ABSOLUTE DEVIATION VARIABLES: 
MO - MEAN DEVIATION 'OVERAGE' 
MU - MEAN DEVIATION 'UNDERAGE' 
MD - MEAN 'TOTAL DEVIATION 
KVO - NUMBER OF 'OVERAGE' ESTIMATES IN THE 
AGGREGATE 
KVU - NUMBER OF 'UNDERAGE' ESTIMATES IN THE 
AGGREGATE 
KVD - NUMBER OF 'TOTAL' ESTIMATION DEVIATIONS 


IN THE 


RR' REPRESENTS A RISK VARIABLE. AN 


RRML (MLE) 
RRJ 


RRTEBM 


AGGREGATE 


THESE PREFIXES ARE APPENDED TO THE BASIC VARIABLES TO 


+ + + + T T+ T+ cb ob cb Gb Gb cb cb cb Gb cb cb cb cb Gb cb Gb Cb cb cb Gb OR Gb Cb Cb cb cb cb cb cb Gb b Gb cb cb cb ck cb cb cb cb Gb πε δ δα πα ο % 


МАК00500 
MAROO510 
МАК00520 
МАК00530 


* MAROO540 


МАК00550 
МАК00560 
МАК00570 
МАК00580 
МАК00590 
MAROO600 
МАК00610 
MARO0620 
МАК00630 
MARO0640 
МАК00650 
МАК00660 
МАК00670 
MAROO680 
МАК00690 
MAROO700 
МАК00710 
МАК00720 
МАК00730 
MARO0 740 
MAR00750 
MAR00760 
MAR00770 
МАК00780 
МАК00790 
МАК00800 
МАК00810. 
МАК00820 
МАК00830 
MAROO840 
MAROO850 
MAROO860 
МАК00870 
МАК00880 
МАК00890 
МАК00900 
МАК00910 
МАК00920 
МАК00930 
МАК00940 
МАК00950 
МАК00960 
МАК00970 
МАК00980 
МАК00990 
МАК01000 
МАК01010 
МАК01020 
МАК01030 
MARO1040 
MARO1050 


+ + + + % + + +f + + + + + + + + 


FORM THE MAD VARIABLES AS FOLLOWS 
MOTSCA MUTSCA MDTSCA 


MOTSAM MUTSAM 
MOPMLE MUPMLE 


HOP J 


MOP JM 
ΠΟΡΟΙ 
MORJTIMIMURTIM 


MUEJ 


MUP IM 
МОРІ 


MDTSAM 
MDPMLE 
MDPJ 
ΠΡΙ 
ΠΠ; 
MDE TIT 


KVOTSA 
KVOTSM 
KVOPML 
KVOPJ 

KVOPJM 
МҮШЕСІ 
KVOP1M 


KVUTSA 
KVUTSM 
KVUPML 
KVUPJ 

KVUPJM 
KVUPJ1 
KVUP1M 


KVDTSA 
KVDTSM 
KVDPML 
KVDPJ 

KVDPJM 
KVDPJ1 
KVDP1M 


(TSCA) 


( TSCAM) 


ΓΕ ΗΕ) 
(PJ) 

(РО) 
(BTT) 
(PJ1M) 


KVDPEB 
KVDPBM 
KVDPTB 
KVDBTM 


KVOPEB 
KVOPBM KVUPBM 
KVOPTB KVUPTB 
KVOBTM KVUBTM 


МОРЕВ МОРЕВ KVUPEB 
МОРЕВМ МОРЕВМ 
IOP TEB HUPTEB 


MOPTBM MUPTBM 


MDPEB 

MDPEBM 
MDPTEB 
MDPTBM 


INPUT/OUTPUT: 


INPUT IS TAKEN FROM A FILE CREATED FROM A TAPE 
PROVIDED BY NPRDC 
OUTPUT IS TO THE PRINTER 


(PEB) 

GEE BM) 
(РТЕВ) 
CELE BM) 


X X X + X c Gb Gb ob ex 


* 


ж 
ж 
ж 
ж 
Үс 


λαο οσους 


FIXED PARAMETERS 
PARAMETER (MXYR=7, MXLOS=31, MXGRD=10, MXMOS=15) 
INPUT PARAMETERS 


жж 


Je 


XX 


XX 


ÚK IPGER ST1, 512, 


LYR 


INTEGER 51051, 51052 

INTEGER SMOS(MXMOS), NMOS 

INTEGER SGRD(MXGRD), NGRD 
VARIABLES IN ORIGINAL SPACE 

INTEGER T, FLAG, IVYR(MXYR) 

INTEGER*2  D(MXMOS, MXLOS, MXGRD) 

INTEGER KV(MXYR) 


REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 


AVCINV,AINV,ATT,ATT1 

CINV(MXMOS , MXLOS ,MXGRD ,MXYR) , YCMXMOS , MXLOS , MXGRD , MXYR) 
ACINV(MXMOS , MXLOS , MXGRD) 

PJ(MXMOS,MXLOS,MXGRD), PJM(MXMOS,MXLOS , MXGRD) 
PMLE(MXMOS,MXLOS,MXGRD), PJLT(MXMOS,MXLOS,MXGRD) 
PJ1(MXMOS , MXLOS ,MXGRD), PJ1M(MXMOS ,MXLOS , MXGRD) 
PEB(MXMOS , MXLOS ,MXGRD), PEBM(MXMOS , MXLOS , MXGRD) 
PTEB(MXMOS , MXLOS ,MXGRD), PTEBM(MXMOS , MXLOS , MXGRD) 
RRTSCA(MXYR),RRIC(MXYR),RRILT(MXYR) , RRML(MXYR) 
В Л, КЕЛЛИ ИХ УВ), ВЕТ СМХУВ), ВВЫСЬ) 
ВЕРЕВОК ТВ ), КВРЕВЫСМХУК) 
RRPTEB(MXYR) , RRPTBM(MXYR) 


VARIABLES IN TRANSFORMED SPACE 


REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 


X(MXMOS , MXLOS , MXGRD , MXYR) 
XT(MXMOS ,MXLOS , MXGRD , MXYR) 
XB(MXMOS,MXLOS,MXGRD), TSCACMXMOS ,MXLOS , MXGRD) 
XTB( MXMOS , MXLOS , MXGRD) 

XJ(MXMOS,MXLOS,MXGRD), XJLT(MXMOS,MXLOS , MXGRD) 
XEB(MXMOS,MXLOS,MXGRD), XJ1(MXMOS,MXLOS,MXGRD) 
XTEB(MXMOS , MXLOS , MXGRD) 

XMLE( MXMOS , MXLOS , MXGRD) 
SHREB(MXMOS , MXLOS , MXGRD) 
SHRTEB( MXMOS , MXLOS , MXGRD) 

RML(MXYR), RSL(MXYR),RXEB(MXYR) 

RXTEB(MXYR) 

RTSCACMXYR) , RTJ( MXYR) ,RTJLT(MXYR) 
RJ(MXYR),RJLT(MXYR), RJI(MXYR) 
TSCAM(MXMOS , MXLOS , MXGRD) 


. 9 
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- 


МАК01060 
MARO1070 
MARO1080 
MARO1090 
MARO1100 
MARO1110 
MARO1120 
МАК01130 
MARO1140 
MARO1150 
MARO1160 
МАК01170 
МАК01180 
МАК01190 
МАК01200 
МАК01210 
MARO1220 
МАК01230 
MARO1240 
MARO1250 
MARO1260 
MARO1270 
MARO1280 
MARO1290 
МАК01300 
MARO1310 
МАКО1320 
MARO1330 
МАК01340 
МАК01350 
МАК01360 
MARO1370 
МАК01380 
МАК01390 
MARO1400 
МАК01210 
MARO1420 
MARO1430 
MARO1440 
МАК01450 
МАК01460 
МАК01270 
MARO1480 
MARO1490 
МАК01500 
MARO1510 
МАК01520 
МАК01530 
МАК01540 
MARO1550 
МАК01560 
MARO1570 
MARO1580 
МАК01590 
МАК01600 
МАК01610 


REAL V(MXMOS ‚MXLOS ‚MXGRD) 
REAL VT(MXMOS ,MXLOS , MXGRD) 
REAL*8 SHR, SHRX, SHR1, SSB, SST, SSE, XBB 
үч VARIABLES TOR FIGURE OF MERIT 
INTEGER KVOTSA(MXYR), KVOTSM(MXYR) , KVOPJ(MXYR), KVOPJ1(MXYR) 
INTEGER KVOPJM(MXYR) ,KVOPEB(MXYR) ,KVOPIN(MXYR) , KVOPBM( MXYR) 
INTEGER KVOPTB(MXYR) ,KVOBTM( MXYR) 
INTEGER KVOPML(MXYR) 
INTEGER KVUTSA( MXYR) , KVUTSM( MXYR) , KVUPJ( MXYR) ,KVUPJ1(MXYR) 
INTEGER KVUPJM( MXYR) , KVUPEB( MXYR) , KVUP1M( MXYR) , KVUPBM( MXYR) 
INTEGER KVUPTB(MXYR) , KVUBTM( MXYR) 
INTEGER KVUPML(MXYR) 
INTEGER KVDTSA(MXYR) , KVDTSM( MXYR) , KVDPJ( MXYR) ,KVDPJ1( MXYR) 
INTEGER KVDPJM(MXYR) ,KVDPEB( MXYR) , KVDP 1M( MXYR) , KVDPBM( MXYR) 
INTEGER KVDPTB(MXYR) , KVDBTM( MXYR) 
INTEGER KVDPML(MXYR) 
REAL MUTSCA(MXYR) , MUTSAM(MXYR) , MUPJ( MXYR) , MUPJ1( MXYR) , MUPJM( MXYR) 
REAL MUPJIM(MXYR) , MUPEBM(MXYR) ,MUPEB( MXYR) , MUPMLE( MXYR) 
REAL MUPTBM(MXYR) , MUPTEB( MXYR) 
REAL MOTSCA(MXYR) , MOTSAM(MXYR) , MOPJ( MXYR) ,MOPJ1( MXYR) , MOPJM( MXYR) 
REAL MOPJIM(MXYR),MOPEBM( MXYR) , MOPEB( MXYR) , MOPMLE( MXYR) 
REAL MOPTBM(MXYR),MOPTEB( MXYR) 
REAL MDTSCA(MXYR) , MDTSAM( MXYR) , MDPJ( MXYR) ,MDPJ1(MXYR) , MDPJM( MXYR) 
REAL MDPJ1M(MXYR),MDPEBM( MXYR) , MDPEB(MXYR) , MDPMLE(MXYR) 
REAL MDPTBM(MXYR),MDPTEB( MXYR) 
REAL MILLER 
** INPUT DATA AREAS 
INTEGER MOS, LOS, GRADE, YR, INV 
INTEGER MOS1,LOS1,GRADE1,YR1,INV1 
XX INPUT INITIALIZATION 
DATA KV/MXYR*0/ 
DATA MOPMLE/MXYR*O / , MUPMLE /MXYR*0/ , MDPMLE /MXYR*0/ 
DATA KVOPML/MXYR*0/ , KVUPML/MXYR*0/ ,KVDPML/MXYR*0/ 
DATA MOTSCA/MXYR*0/ ,MUTSCA/MXYR*0/ ,MDTSCA/MXYR*0/ 
DATA KVOTSA/MXYR*0/ ,KVUTSA/MXYR*0/ ,KVDTSA/MXYR*0/ 
DATA MOTSAM/MXYR*0/ , MUTSAM/MXYR*0/ , MDTSAM/MXYR*0/ 
DATA KVOTSM/MXYR*0/ , KVUTSM/ MXYR*0/ , KVDTSM/MXYR*0/ 
DATA MOPJ/MXYR*0/ ,MUPJ/MXYR*0/ ,MDPJ/MXYR*0/ 
DATA KVOPJ/MXYR*0/ ,KVUPJ/MXYR*0/ ,KVDPJ/MXYR*0/ 
DATA MOPJM/MXYR*0/ ,MUPJM/MXYR*0/ ,MDPJM/MXYR*0/ 
DATA KVOPJM/MXYR*0/ , KVUP JM/MXYR*O / ,KVDPJM/MXYR*0/ 
DATA MOPJ1/MXYR*0/ ,MUPJ1/MXYR*0/ ,MDPJ1/MXYR*0/ 
DATA KVOPJ1/MXYR*0/ ,KVUPJ1/MXYR*0/ , KVDPJ1/MXYR*0/ 
DATA MOPJ1M/MXYR*0/ ,MUPJ1M/MXYR*0/ ,MDPJ1M/MXYR*0/ 
DATA KVOP1M/MXYR*0/ , KVUP1M/MXYR*0/ , KVDP1M/MXYR*0/ 
DATA MOPEB/MXYR*0/ ,MUPEB/MXYR*0 / ,MDPEB/MXYR*0/ 
DATA KVOPEB/MXYR*0/ ,KVUPEB/MXYR*0/ ,KVDPEB/MXYR*0/ 
DATA MOPEBM/MXYR*0/ ,MUPEBM/MXYR*0/ ,MDPEBM/MXYR*0/ 
DATA KVOPBM/MXYR*0/ , KVUPBM/MXYR*0/ , KVDPBM/MXYR*0/ 
DATA MOPTEB/MXYR*0/ ,MUPTEB/MXYR*0/ ,MDPTEB/MXYR*0/ 
DATA KVOPTB/MXYR*0/ ,KVUPTB/MXYR*0/ ,KVDPTB/MXYR*0/ 
DATA MOPTBM/MXYR*0/,MUPTBM/MXYR*0/ ,MDPTBM/MXYR*0/ 


DATA KVOBTM/MXYR*0/ ,KVUBTM/MXYR*0/ , KVDBTM/MXYR*0/ 
* 


ποσο ο ο ο ος κ ο ο ος σος κο λος κος κο ο ο ος ος οκ ο ο κκ κκ ος ος Κε ης 


== Ν 
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MARO1620 
MARO1630 
MARO1640 
MARO1650 


` MAR01660 


MAR01670 
MAR01680 
MAR01690 
MAR01700 
MAR01710 
MAR01720 
MAR01730 
MARO1740 
MARO1750 
MARO1760 
МАК01770 
МАК01780 
МАК01790 
MARO1800 
MARO1810 
MARO1820 
MARO1830 
МАК01840 
МАК01850 
MARO1860 
МАК01870 
МАК01880 
МАК01890 
MARO1900 
MARO1910 
MARO1920 
MARO1930 
MARO1940 
MARO1950 
MARO1960 
МАК01970 
МАК01980 
МАК01990 
МАК02000 
МАК02010 
MARO2020 
MARO2030 
MARO2040 
MARO2050 
MARO2060 
MARO2070 
MARO2080 
MARO2090 
MARO2100 
MARO2110 
MARO2120 
MARO2130 
MARO2140 
MARO2150 
MARO2160 


XX READ PARAMETER DATA TO USE IN THIS RUN. ж 


ET Tee re 
МВТТЕ( 5,100) 
100 FORMAT(' ENTER 1ST YEAR AND LAST YEAR TO USE FOR ESTIMATION') 
READ(5,*) ST1,ST2 
WRITE(6,*) 'ESTIMATION YEARS: ',ST1,ST2 
ж 
WRITE(5,101) 
101  FORMAT(' ENTER LAST YEAR PRESENT IN THE DATA BASE') 
READ(5,*) LYR 
WRITE(6,*) 'LAST YEAR AVAILABLE FOR VALIDATION: ',LYR 
ж 
МЕІТЕ(5,102) 
102  FORMAT(' ENTER NO. OF MOS FOLLOWED BY ARRAY OF SAME LENGTH') 
READ(S,*) NMOS, (SMOS(I), I21,NMOS) 
WRITE(6,*) 'MOS SELECTED: ', (SMOS(I), I=1,NMOS) 
* 
МЕІТЕ(5,103) 
103  FORMAT( ENTER 1ST AND LAST LOS VALUE TO USE') 
READ(5,*) SLOS1,SLOS2 
WRITE(6,*) 'LOS RANGE: * ,SLOS1,SLOS2 
ж 
WRITE(5,104) 
104 . FORMAT(' ENTER NO. OF GRADE FOLLOWED BY ARRAY OF SAME LENGTH') 
READ(5,*) NGRD, (SGRD(I), I=1,NGRD) 
WRITE(6,*) 'GRADES SELECTED', (SGRD(1), I=1,NGRD) 
ж 
МЕІТЕ(5,105) 
105 ЕОЕМАТ( ENTER DEE FACTOR’) 
READ(5,*) DEE 
WRITE(6,*) 'DEE FACTOR: ' ,DEE 
* 
жж COMPUTE START AND END YEARS FOR VALIDATION 
NYR=LYR-ST1+1 
NVYR=LYR -ST2 
NLOS=SLOS2-SLOS1+1 
T=ST2-ST1+1 


DO 41 I=1,NVYR 
IVYR(1)=ST2+1 
41 CONTINUE 
ee 


= INITIALIZE D ARRAY TO ALLOW FINDING MISSING VALUES. XX 
er 
DO 1 I=1,MXMOS 
DO 2 J=1,MXLOS 
DO 3 K=1,MXGRD 


D(1,J,K)=-9999 

3 CONTINUE 
2 CONTINUE 
1 CONTINUE 
RITTER ra ee 
** ` TAPE PROCESSING TO INPUT INVENTORY AND ATTRITION VALUES, ** 
XX DEVELOPE D MATRIX ACCORDING TO EXTRACTION CRITERIA. ** 
πο οσο ος κακο κοκ οσο ο οκ οκ е а а а ааа 

ҒІАС-0 

IEOF=0 


ა.დ 
` 
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MAR02170 
MAR02180 
МАК02190 
MAROZINN 
ШАҚ 2210 
MARO2220 
MARO2230 
MARO2240 
МАКО2250 
MARO2260 
MARO2270 
MARO2280 
MARO2290 
MARO2300 
МАК02310 
MARO2320 
МАКО2330 
MARO2340 
МАК02350 
МАК02360 
МАК02370 
MARO2380 
МАК02390 
MARO2400 
MARO2410 
MARO2420 
МАКО2430 
МАК02440 
MARO2450 
MARO2460 
МАК02%70 
MARO2480 
MARO2490 
MARO2500 
MARO2510 
MARO2520 
MARO2530 
MARO2540 
MARO2550 
MARO2560 
MARO25 70 
MARO2580 
MARO2590 
MARO2600 
MARO2610 
МАК02620 
MARO2630 
MARO2640 
MARO2650 
MARO2660 
MARO2670 
MARO2680 
MARO2690 
MARO2700 
МАК02710 
MARO2720 : 


сус 


12 


ro 


ee READER(MOS1,L0S1,GRADE1,YR1,INV1,ATT1, NMOS,SMOS, 
NGRD,SGRD, 51051, 51052, u ST2 IM1, TEMIGI, ГТ ТЕОЕ) 


ТЕ(ТЕОЕ . NE. O) THEN 
WRITE(6,*) '**®** NO DATA MEETS SELECTIONS REQS' 
STOP 

ENDIF 


IF(IEOF . NE. Ὁ) GO TO 9 

MOS=MOS1 

LOS=LOS1 

GRADE=GRADE1 

YR=YR1 

ІМУ-ІМУІ 

АТТ-АТТ1 

IM=IM1 

IL-IL1 

IG=IG1 

IT=IT1 

IF(YR. GT. ST1) FLAG - 

CHECK CASE WHERE DATA BEGINS PAST 15Т YR 

IF(YR .GT. ST1) CINV(IM,IL,IG,IT-1)=. 5*FLOAT( INV) 
CALL READER(MOS1,LOS1,GRADE1,YR1,INV1,ATT1, NMOS,SMOS, 


uj NGRD,SGRD, SLOS1,SL105277 5718272, IM μι Ποπ. ο» 


ONLY FOR CHECKING PURPOSES. MARK D TO INDICATE SOME YR PRESENT 
D(IM,IL,IG)=0 


IF(.NOT. (IEOF.EQ.O .AND. MOSI1. EQ. MOS . AND. LOS1. EQ. LOS 
. AND. GRADE1. EQ. GRADE)) GO TO 8 
CENTRAL INV. FOR YR1 - 1 
III II) ცი...) I.) 
CINV(1M1,1L1,1G1,IT1-1)=. 5*FLOAT(INV1) 
AINV=. 5*FLOAT( INV) 
I=YR+1 
FLAG=1 
ELSE 
AINV=. 5*FLOAT( INV+INV1) 
END IF 
CINV(1M,IL,IG,IT)=AMAX1(AINV, ATT) 
YCIM,IL,IG,IT)-ATT 
YR=YR1 
INV=INV1 
ΑΤΤΞΑΤΤΙ 
ІМ-ІМІ 
IL=IL1 
IG=IG1 
ІТ-ІТІ 
CALL READER(MOS1,LOS1,GRADE1,YR1,INV1,ATT1, NMOS,SMOS, 


* NGRD,SGRD, SLOST,SLOS27 ST1E ST2 oi TIS TEGEM 


GO TO 12 
CONTINUE 
WHEN YEARS MISSING AT THE END и. « 
ТСК. т. სს) πὴ 
FLAG=1 `` 
AINV=. S*FLOAT( INV) 
ELSE 


v? 
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MARO2730 
MARO2740 
MARO2750 
MARO2760 
MARO2770 
MARO2780 
МАК02790 
MARO2800 
МАК02810 
МАК02820 
МАК02830 
MARO2840 
MARO2850 
MAR02860 
MAR02870 
MAR02880 
MAR02890 
MAR02900 
MAR02910 
MAR02920 
MAR02930 
MARO2940 
MARO2950 
MARO2960 
МАК02970 
MARO2980 
МАК02990 
MARO3000 
MARO3010 
MARO3020 
MARO3030 
MARO3040 
MARO3050 
MARO3060 
MARO3070 
MARO3080 
MARO3090 
МАК03100 
MARO3110 
MARO3120 
MARO3130 
MARO3140 
МАК03150 
МАК03160 
MARO3170 
MARO3180 
МАК03190 
МАК03200 
MARO3210 
MARO3220 
MARO3230 
MARO3240 
MARO3250 
MARO3260 
МАК03270 
MARO3280 


AINV=FLOAT( INV) 
ED IF 
CINV(IM,IL,IG,IT)=AMAX1(AINV,ATT) 
ШЕНІ, ІТ,,ІС,ІТ)-АТТ 


ВУ ТМ, ТЪ,ТС,ТТ) 
ВНЕ. +2. ХАТТ/С1. +СС) 

' τ᾽. 2.*(1. +АТТ)/(1. +СС) 
CC=ABS( TEMP) 

CC1-ABS(TEMP1) 


ШИЕ ОЛ. 1. .OR. CCl .GT. 1.) WRITE(6,*) '*** TEMP, TEMP1=', 
Ж TEMP,TEMP1, '***MOS,LOS,GR,YR=' ,SMOS(IM),IL,SGRD(IG) ,IT 
@Ш ТО 11 
CONTINUE 
AAA ο λος ΕΣΣΟ κλπ κκ ο ος λος πο ΚΚ 
κ CALCULATE AVERAGE CENTRAL INVENTORY AND MAXIMUM XX 
к LIKELIHOOD ESTIMATE. Je 


AAA RRA AAA AA 


DO 10 IM=1,NMOS 
DO 20 IL=1,NLOS 
DO 30 IG=1,NGRD 
SUMY=0 
ACINV( IM, IL, 1G)=0 
"„> 35 II=1 L 
SUMY=SUMY+Y(IM,IL,IG,IT) 
ESENVCIM,ID,IG)-ACINV(CIM,IL, IG)TCINV(CIM, IL, IG, IT) 
5 CONTINUE 
ШҒАСТЫУ( ІМ, ІІ,,ІС) .МЕ. 0) THEN 
PMLE( IM, IL, 1G)=SUMY/ACINV( IM, IL, IG) 
EESE 
PMLE( IM, IL, 1G)=0 
ENDIF 
ACINV(IM, IL, 1G)=ACINV(1M, IL,IG)/T 
DO 40 IT=1,NYR 


x IF ANY CENTRAL INVENTORY>0 (OVER EST. YRS T) THEN D=1 
ФҮСІНУСІМ, ІІ, ІС,ІТ). 6І.0 .AND. IT. LE.T) DCIN,IL,IG)=1 

r IF ANY CENTRAL INVENTORY>O (OVER VAL. YRS > T) THEN COUNT 
КОО (IH I IL IG IT).GT.0 . AND. IT.GT.T) KV(IT=T)=KV( IT-T)+1 


ü TRANSFORMATION OF CINV USING THE FREEMAN-TUKEY 

DOUBLE ARCSIN FORMULA. 
ENDS DOODUCIMOIDSIGSIT)SYCIM,ID,IG, IT), XCIM, IL, IG,IT)) 
BEN CT IGSIT)EXCIM;IB,IG,IT) / SQRICO. STCINV(CIM,IL, IG, IT)) 


ж 


0 CONTINUE 

0 CONTINUE 

0 CONTINUE 

0 CONTINUE 

TÜTE TE TE TA TE TE Te ee ee Te de Te УГУС 

к CHECK FOR MISSING DATA, COMPUTE THE TRANSFORMED XX 

к CELL AVERAGES AND GRAND MEAN. XX 

АКК еее еее еее СТК КАКАСИ СУС Ж 
ХВВ-0 : 


KK=0 

DO 50 IM=1,NMOS С 

DO 60 IL=1,NLOS 
DO 70 IG=1,NGRD 


y 
4.” 
` 
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ЧАК03290 
МАК03300 
ЧАК03310 
MARO3320 
MARO3330 
MARO3340 
MARO3350 
MARO3360 
MARO3370 
MARO3380 
MARO3390 
MARO3400 
MARO3410 
MARO3420- 
MARO3430 
MARO3440 
MARO3450 
MARO3460 
MARO3470 
MARO3480 
MARO3490 
MARO3500 
MARO3510 
MARO3520 
MARO3530 
MARO3540 
МАК03550 
МАК03560 
МАК03570 
MARO3580 
MARO3590 
MARO3600 
MARO3610 
MARO3620 
MARO3630 
MARO3640 
MARO3650 
МАК03660 
МАК03670 
МАК03680 
МАК03690 
MARO3700 
MARO37 10 
МАК03720 
MARO3730 
MARO3740 
МАК03750 
МАК03760 
MARO3770 
MARO3780 
МАК03790 
MARO3800 
MARO3810 
MARO3820 
МАК03830 
MARO3840 


XX 


80 


XX 


ж 


70 
60 
50 


rer ee TAT AT A I I TONE TEN IE 


BAYES VARIANCE ITERATION SECTION. 
nee 


er 


XX 
#99 


522 


220 
510 
500 


530 


393 


532 
531 


ІРСРСІМ,ТЫ, 16) . ЕО. THEN 


REPORT MISSING COMBINATION IM,IL,IG & CLEAR D ENTRY 


END IF 


D(IM,IL,IG)=0 
ILD=SLOS1+IL-1 
FLAG=1 


KK=KK+D(1M, IL, 1G) 


ВИМ, 


IL,IG)=0 


XTB(IM,IL,IG)=0 


DO 80 


XBCIM,IL, IG) SABCIN, IL IG EXC IN IL TOIT) 
XTB(IM,IL,1G)=XTBCIM, TL IG) ИИСИ ШЕН 


IT=1,T 


CONTINUE 


ΠΠ. 


IL,IG)=XB(IM,IL,IG)/T 


XTB( IM .,IL,.IG)=XTB( IM ID IC) YI 


EMPIRICAL BAYES PREPARATION ON THE TRANSFORMED SCALE. 


XEB(IM,IL,IG)=XB(IM,IL,IG) 
XTEB( IM,IL,IG)-XTB(IM,IL, IG) 


XBB=XBB+D( IM, IL, 1G)*XB(1M, IL, 1G) 


CONTINUE 


CONTINUE 


CONTINUE 


XBB=XBB/KK 


ІС-0 


COMPUTE TRANSFORMED SCALE CELL VARIANCE. 


DO 500 IM-1,NMOS 


DO 510 


IL=1,NLOS 


DO 520 IG=1,NGRD 
IF (D( IM, ID IG) ONE ТІ) СОТО#520 
AVSQR = 0.0 


DO 


522 IT=1,T 


AVSQR = AVSQR + SQRT(0.5 + CINV(IM,IL,IG,IT)) 
Z = XEB(IM,IL,IG) + (AVSQR/T) * 1.5708 


V(IM,IL,IG) = VAR(Z) 
CONTINUE 
CONTINUE 


CONTINUE 


A=0 
AI=A 
SALPH=0 


DO 531 IM=1,NMOS 


DO 532 


IL=1,NLOS 


DO 533 IG=1,NGRD 


IF (D(IM,IL,IG). EQ. 1) SALBH 2 SALEN TET A ПС ШЕТІ) 


CONTINUE 


CONTINUE 


28=0 


COMPUTE BAYES ALPHA AND GAMMA PARAMETERS 


DO 537 IM=1,NMOS 


DO 538 


IL-1,NLOS 


a 
` + 
` 


58 


MAR03850 
MAR03860 
MAR03870 
MARO3880 
MARO3890 
MARO3900 
MARO3910 
MARO3920 
MARO3930 
MARO3940 
MARO3950 
MARO3960 
MARO3970 
MARO3980 ` 
MARO3990 
MARO4000 
МАК04010 
МАК04020 
MARO4030 
MARO4040 
MARO4050 
MARO4060 
МАК04070 
МАК04080 
MARO4090 
MARO4100 
MARO4110 
МАК04120 
МАК04130 
ΜΑΚΟΦΙ40 
МАКО4 150 
MARO4160 
MARO4170 
MARO4180 
MARO4 190 
MARO4 200 
MARO4210 
МАК04220 
МАК04230 
МАК04240 
MARO4250 
MARO4260 
MARO4270 
MARO4280 
МАК04290 
MARO4300 
МАК08310 
MARO4320 
MARO4330 
MARO434O 
MARO4350 
MARO4360 
MARO4 370 
MARO4 380 
MARO4 390 
МАК04400 


539 
538 
637 


542 
541 


540 
XX 


550 


590 
580 
570 


599 


Je 


602 


DO 539 IG=1,NGRD 
IF (D(IM,IL,IG).EQ. 1) THEN 
ALPH 2 1./(A + V(IM,IL,IG)) 
GAM = ALPH / SALPH 
28 ZB + GAM*XB(IM,IL,IG) 
ENDIF 
CONTINUE 
CONTINUE 
CONTINUE 
F=0 
G=0 
DO 540 IM=1,NMOS 
DO 541 IL=1,NLOS 
DO 542 IG=1,NGRD 
IF (D(IM,IL,IG).EQ. 1) THEN 
ALPH = 1./(A + V(IM,IL,IG)) 


5 = F + (ALPH * (XB(IM,IL,IG) - ZB)**2) 
G = G + (ALPH * (XBCIM,IL,IG)-ZB))**2 
ENDIF 
CONTINUE 
CONTINUE 
CONTINUE 


COMPUTE THE PRIOR VARIANCE (A) 
À = À - (KK-1-F)/G 


IF (A.LE.0.) THEN 
A=0 
GOTO 550 
ENDIF 
IF (ABS(A-AI).GT..0001) GOTO 530 
COMPUTE THE EMPIRICAL BAYES ATTRITION RATES 
DO 570 IM=1,NMOS 
DO 580 IL-1,NLOS 
DO 590 IG=1,NGRD 
IF (D(IM,IL,IG).EQ. 1) THEN 
DEN = A + V(IM,IL,IG) 
XEB(IM,IL,IG) = (A/DEN)*XB(IM,IL,IG)+((DEN-A)/DEN )*ZB 
ENDIF 
CONTINUE 
CONTINUE 
CONTINUE 
IC=IC+1 
IF (A.EQ.0) GOTO 599 
ШІ (ІС. ІТ. 10) GOTO 499 


CONTINUE 


COMPUTE THE EMPIRICAL BAYES SHRINKAGE RATES 
DO 600 IM=1,NMOS 
DO 601 IL=1,NLOS 

DO 602 IG=1,NGRD | 
SHREB(IM,IL,IG) = O | : 
IE DO TI IL, ,IG). EQ. 1) THEN 

SHREB(IM,IL,IG) = V(IM,IL,IG)/(A + V(IM,IL,IG)) 
ENDIF 
CONTINUE 


>. 


59 


MARO4410 
MARO4420 
MARO4430 
МАК04%50 
MARO4450 
МАК04460 
МАК04470 
MARO4480 
MARO4490 
МАК04500 
МАК04510 
MARO4520 
MARO4530 
MARO4540 | 
MARO4550 
MARO4560 
MARO4570 
MARO45 80 
МАК04590 
МАК04600 
МАК04610 
MARO4620 
MARO4630 
MARO4640 
MARO4650 
МАК04660 
MARO4670 
MARO4680 
MARO4690 
MARO4700 
МАК04710 
MARO4 720 
MARO4 730 
MARO4740 
MARO4750 
MARO4760 
MARO4770 
MARO4780 
MARO4 790 
МАК04800 
MARO4810 
MARO4820 
MARO4830 
MARO4840 
MARO4850 
MARO4860 
МАК08870 
MARO4880 
MARO4890 
MARO4900 
МАК04910 
МАК04920 
МАК08930 
МАК04940 
МАК04950 
МАК04960 


601 
800 


эзе 


Зе СКК КК еее еее HEE Te FETE TE AA AAA 


XX 
699 


723 


720 
710 
700 


730 


733 
792 
731 


739 
738 
737 


742 
781 


CONTINUE 


CONTINUE 
ee 


ALTERNATE BAYES VARIANCE ITERATION SECTION. 
COMPUTE TRANSFORMED SCALE CELL VARIANCE. 


DO 700 IM=1,NMOS 


DO 710 IL=1,NLOS 
DO 720 IG=1,NGRD 
VT(IM,IL,IG) = 0.0 
IF (D(IM,IL,IG).NE. 1) GOTO 720 
DO 725 IT=1,T 
ZT = X(IM,IL,IG,IT) + 1.5708 * SQRT(0.5 + 
CINV(IM,IL,IG,IT)) 
VI(IM,IL,IG) s VI(IM,IL,IG) * VAR(ZT)/(0.5 - 
CINV(IM,IL,IG,IT)) 
VT(IM,IL,IG) = УТІМ ID. IG) Fo 
CONTINUE 
CONTINUE 
CONTINUE 


AT=0 
ATI=AT 
SALPHT=0 
DO 731 IM=1,NMOS 
DO 732 IL=1,NLOS 
DO 733 IG=1,NGRD 


IF (D(IM,IL,IG).EQ. 1) SALPHT= SALPHT+ 1./( AT ΥΠΗΡ 9 


CONTINUE 
CONTINUE 
ZTB=0 
COMPUTE BAYES ALPHA AND GAMMA PARAMETERS 
DO 737 IM=1,NMOS 
DO 738 IL=1,NLOS 
DO 739 IG=1,NGRD 
IF (D(IM,IL,IG).EQ.1) THEN 
ALPHT= 1. /(AT+ VT(IM,IL,IG)) 
GAMT = ALPHT / SALPHT 
ZTB - ZTB* GAMT*XTB( IM,IL,IG) 
ENDIF 
CONTINUE 
CONTINUE 
CONTINUE 
PZTB=0. 5*(1. *SIN(ZTB)) 
FT=0 
GT=0 
DO 740 IM=1,NMOS 
DO 741 IL=1,NLOS 
DO 742 IG=1,NGRD 
IF (D(IM,IL,IG).EQ. 1) THEN 
ALPHT= 1. /(AT * VI(IM,IL,IG)) 


FT = FT + (ALPHT * (XTB(IM,IL,IG) = 2:8)XX2) 
СІ = GT + (ALPHT * (XTB(IM, IL, 16) е 
ENDIF წ. 
CONTINUE 
CONTINUE 


` 
~P 
› 


60 


МАКО5970 
МАК04980 
МАК04990 
МАК05000 
MARO5010 
МАК05020 
МАК05030 
MAROS049 
MARO5050 
МАК05060 
МАКО5070 
МАКО5080 
MAR05090 
MARO5100' 
MARO5110 
МАК05120 
МАК05130 
МАК05140 
МАК05150 
МАКО5160 
МАКО5170 
МАКО5180 
МАКО5190 
МАКО5200 
МАКО5210 
МАК05220 
МАК05230 
МАК05240 
МАКО5250 
МАКО5 260 
МАКО5 270 
МАКО5280 
MARO5 290 
MARO5 300 
MARO5 310 
MARO5320 
MARO5330 
MARO5340 
МАК05350 
MARO5360 
МАК05370 
МАК05380 
МАК05390 
МАК05400 
MARO5410 
МАК05420 
МАК05430 
МАК05440 
МАК05450 
МАК05460 
МАК05470 
МАК05480 
МАК05490 
МАК05500 
МАК05510 
MARO5520 


740 
жж 


жж 
750 


790 
780 
770 


802 
801 
800 


CONTINUE 
COMPUTE THE PRIOR VARIANCE (A) 
AT= AT- (KK-1-FT)/GT 


i) (AT. LE. 0, ) THEN 
AT=0 
GOTO 750 
ENDIF 
ШЕ (АВ5(АТ-АТІ). СТ.. 0001) СОТО 730 
COMPUTE THE EMPIRICAL BAYES ATTRITION RATES 
DO 770 IM=1,NMOS 
DO 780 IL=1,NLOS 
DO 790 IG=1,NGRD 
IF (D(IM,IL,IG).EQ. 1) THEN 
DENT = AT + VI(IM,IL,IG) 
XTEB(IM,IL,IG) = (AT/DENT)*XTB(IM,IL,IG) + 
( (DENT-AT) /DENT)*ZTB 
ENDIF 
CONTINUE 
CONTINUE 
CONTINUE 
COMPUTE THE EMPIRICAL BAYES SHRINKAGE RATES 
DO 800 IM=1,NMOS 
DO 801 IL=1,NLOS 
DO 802 IG=1,NGRD 
SHRTEB(IM,IL,IG) = 0 
IF (D(IM,IL,IG).EQ. 1) THEN 
SHRTEB(IM,IL,IG) = VT(IM,IL,IG)/C(AT + VT(IM,IL,IG)) 
ENDIF 
CONTINUE 
CONTINUE 
CONTINUE 


Tr FA root 


XX 


JAMES-STEIN SHRINKAGE RATES. XX 


AAA a rr 


Kr 


120 


110 
95 
90 


XX 


SST=0 
SSE=0 
SSB=0 
COMPUTE THE SUMS OF SQUARES (SST, SSE AND SSB) 
DO 90 IM=1,NMOS 
DO. 95 IL=1,NLOS 
DO 110 IG=1,NGRD 
TEMP1=0 
TEMP2=0 
DO 120 IT=1,T 
TEMP1=TEMP1+(X( IM, IL, 1G, 1T)-XBB)**2 
TEMP2=TEMP2+(X( IM, IL, IG, IT)-XB( IM, IL, IG) )**2 
CONTINUE 
SST=SST+D( IM, IL, 1IG)*TEMP1 
SSE=SSE+D( IM, IL, IG)*TEMP2 
CONTINUE 
CONTINUE : 
CONTINUE 
SSB-SST-SSE 
COMPUTE THE JAMES-STEIN SHRINKAGE RATES 
SHRX=(SSE/SSB)* (KK-3)/(KK*(T-1)+2) 


. Ф 


61 


MAR05530 
МАК05540 
МАК05550 
МАК05560 


` HAR05570 


MAR05580 
MAR05590 
MAR05600 
MAR05610 
MARO5620 
МАВ05630 
MARO5640 
MAROS650 
MAROS 660 
MAROS670 
MAROS680 
MAROS690 
MAROS 700 
MARO5710 
ΜΑΒΟ5720 
MARO5730 
MAROS 740 
MARO5750 
MARO5760 
МАВ05770 
MARO5780 
MAROS790 
MARO5800 
MARO5810 
MARO5820 
MARO5830 
MARO5840 
MARO5850 
MARO5860 
MARO5870 
MARO5880 
MARO5890 
MARO5900 
MAROS910 
ΜΑΚΟ5920 
MARO5930 
MAROS940 
MAROS950 
MARO5960 
MARO5970 
ΜΑΒΟ5980. 
MARO5990 
MARO6000 
MARO6010 
MARO6020 
MARO6030 
MARO6040 
MARO6050 
MARO6060 
MARO6070 
MARO6080 


SHR=DMINI(SHRX, 


1. D0) 


SHR=DMAX1(SHR , 0.D0) 


SHR1=T*SHR 


ee 
LIMITED TRANSLATION OPTION & JAMES-STEIN 


XX 
Je 


ON TRANSFORMED SCALE. 


XX 
deve 


Fete eet ere Tove deve AH TENTH HATE fefe efe eee dene e de e fee e de eje defe efe fee fede eee ee ITER IIH WICH 
NN-KK*(T-1) 
DO 125 IM=1,NMOS 
DO 130 IL-1,NLOS 
DO 140 IG-1,NGRD 
XJ(IM,IL,IG)-0 
XJ1(IM,IL,IG)-0 
IF (D(IM,IL,IG).EQ. 1) THEN 
XJC IM, IL, IG)=XBB+( 1-SHR)*(XBC IM, IL, IG) -XBB) 
XJ1( IM, IL, IG)=XBB+(1-SHR1)*(XB(IM,IL,IG)-XBB) 


C w 

ც XXX 
* 

140 


130 
125 


XJLT(IM, IL, 1G)=XBB+(1-SHR*RHO(XB( IM, 1L,IG),XBB,SSE, 
SSB ,DEE ,KK,T,NN))*(XB(IM,IL,IG)-XBB) 
XMLE(IM,IL,IG)=SQRT(.5+ACINV(IM,IL,IG)) 
*ASIN( 2*PMLEC IN, IL,IG)=1) 
ENDIF 
CONTINUE 
CONTINUE 
CONTINUE 


ჯCას წილს 2 729 9 X XV XX XC XX ο ΚΚ ΚΚ κκ κοκ κκ κ Te MCT E е TE ICE FEN FETE FETE 


Ar 


VALIDATION. 


RISK EVALUATED IN TRANSFORM SPACE. 


XX 


HICH PATE AH HE HAE IT NTE HK ITT HICH NIH TTI TT AVE RIT ICT ITE INTE TICK NTE THICK IIH ICE IIIT 
DO 150 IT=1, NVYR 


RTSCA(IT)-20 


RI(IT)=0 
RJ1( IT)=0 
RXEB(IT)=0 
RXTEB( IT)=0 
C week RILT(IT)=0 
RML( IT)=0 
DO 160 IM=1,NMOS 
DO 170 IL=1,NLOS 
DO 180 IG=1,NGRD 
IF(D(IM,IL,IG).EQ.1 . AND. CINV(IM,IL, IG, IT+T). NE. 0) THEN 
RTSCA( IT)=RTSCA(IT)+(X(IM,IL,IG,T+IT)- XB(CIM,IL,IG))**2 
RJ(IT)  RJ(IT(X(IM,IL,IG,T-IT)- XJ(IM,IL,IG))**2 
RJICIT)= RJ1CIT)+(X( 1M, IL,IG,T+IT)- XJ1(1IM,IL,IG))**2 
RXEB( IT)=RXEB( IT)+(X(IM,IL,IG,T+IT)- XEB(IM,IL,IG))**2 
RXTEB( IT)=RXTEB( IT)+((X( IM, IL,IG,T+IT)/SQRT(0.5 + 
ж CINV(IM,IL,IG,IT4T))- XTEB(IM,IL,IG))**2) 
* /VT( IM,IL,IG) 
C жж RJLT (IT)5 RILT(IT)+(X(C IM, IL, IG, T+IT) -XJLT( IM, IL, IG) )**2 
RML (IT)= RML (IT)+(X(IM,IL,IG,T+IT) -XMLE( IM, IL, IG) )**2 
ENDIF 
180 CONTINUE 
170 CONTINUE 
160 CONTINUE | 


RTSCA(IT)=RTSCA(IT)/KV(IT) 
RJ (IT)SRJ (IT)/KV(IT) 
RJ1 (IT)=RJ1 (ІТ)/КУ(ІТ) 
RXEB(IT)=RXEB(IT)/KV(IT) 
C 
|. 62 


` 


MAR06090 
MAR06100 
MAR06110 
MARO6120 
МАК06130 
МАК06140 
МАК06150 
МАК06160 
МАК06170 
МАК06180 
МАК06190 
МАК06200 
МАК06210 
МАК06220 
МАК06230 
MARO6240 
МАК06250 
MARO6260 
MARO6270 
MARO6280 
МАК06290 
МАК06300 
MARO6 310 
МАК06320 
МАК06330 
MARO6340 
МАК06350 
МАК06360 
МАК06370 
МАК06380 
МАК06390 
МАК06400 
MARO6410 
МАК06420 
МАК06430 
MARO6440 
MARO6450 
МАК06460 
МАК06470 
МАК06480 
МАК06490 
МАК06500 
МАК06510 
МАК06520 
МАК06530 
МАК06540 
МАК06550 
МАК06560 
МАК06570 
МАК06580 
МАК06590 
МАК06600 
МАК06610 
MARO6620 
МАК06630 


MARO6640 : 


RXTEB( IT)2RXTEB( IT)/(KVCIT)*T) 
C ***  RILT(IT)=RILT(IT)/KV(IT) 
RML (IT)-2RML (IT)/KV(IT) 
RTJCIT)-RTSCA(IT)/RJCIT) 
С *** RTJLT(IT)-RTSCACIT)/RJLT(IT) 
C *** RSL(IT)-(RJCIT)-RTSCACIT))/CRJLT(CIT)-RTSCA(CIT)) 
150 CONTINUE 
dese ede defe de de ve je ve efe jefe fee e deve dee e ee de ede e eee e ee de efe ve de efe fee dee defe defe e dede delete dede dede dede de dede de 


## DO RISK ANALYSIS ON THE ORIGINAL SCALE. XX 
ee ee Tee Te 
DO 200 IM=1,NMOS 
DO 210 IL=1,NLOS 
DO 220 IG=1,NGRD 
жж AVERAGE CENTRAL INVENTORY OVER ESTIMATION SET 
CINVS=0 
DO 230 IT=1,T 
AVCINV=AVCINV + CINV(IM,IL,IG,IT) 
230 CONTINUE 
AVCINV=AVCINV/T 
TSCA( IM, IL,IG)=D(IM,IL,IG)*SCINV( XB(IM,IL,IG), AVCINV) 
PJ(IM,IL,IG)=D(IM,IL,IG)*SCINV( XJ(IM,IL,IG), AVCINV) 
PJ1( IM, IL,IG)=D(IM,IL,IG)*SCINV( XJ1(IM,IL,IG), AVCINV) 
PEB(IM,IL,IG)=D(IM,IL,IG)*SCINV( XEB(IM,IL,IG), AVCINV) 
PTEB(IM,IL,IG)=D(IM,IL,IG)*0.5*(1. + SIN(XTEB(IM,IL,IG))) 
6 ** — PJLT(IM,IL,IG)SD(IM,IL,IG)*SCINV(XJLT(IM,IL,IG), AVCINV) 
IF (D(IM,IL,IG).EQ. 1) THEN 
TSCAM( IM, IL, IG)=MILLER(  XB(IM,IL,IG), AVCINV ) 

PJM( IM, IL,IG)=MILLER(  XJ(IM,IL,IG), AVCINV ) 
PJIM(IM,IL,IG)=MILLER( XJ1(IM,IL,IG), AVCINV ) 
PEBM(IM,IL,IG)=MILLER( XEB(IM,IL,IG), AVCINV ) 

XMIL = XTEB(IM,IL,IG) * SQRT(0.5 + AVCINV) 
PTEBM(IM,IL,IG)=MILLER( XMIL , AVCINV ) 
ENDIF | 
220 CONTINUE 
210 CONTINUE 
200 CONTINUE 
ee ee ee ee ee re deese 


je VALIDATION. RISK EVALUATED IN ORIGINAL SPACE. Ye 
qyeeseye vede vede eze vec Ve Tete Te ve FOF ὕς Ὃς Ὑς Ὑς усуе УСЭ TOT Te TOTO TOT FE FE TOE IE TC TCT T CTE TC TOV FETE TC FONT ET TORN TOT TTI Te Te I 
DO 240 IT=1, NVYR 
RRTSCA( IT)=0 
RRTSAM( IT)=0 
RRJ( IT)=0 
RRJM( IT)=0 
RRPEB( IT)=0 
RRPTEB( IT)=0 
RRPEBM( IT)=0 
RRPTBM( IT)=0 
КЕЈ1(ІТ)=0 
RRJ1M( IT)=0 
C ** RRILT(IT)=0 | 
RRML( IT)=0 5 | 
KP=0 | 
DO 250 IM=1,NMOS 
DO 260 IL=1,NLOS 
DO 270 IG=1,NGRD 


Ф 
63 


МАК06650 
МАК06660 
МАК06670 
МАК06680 


` MAR06690 


MAR06700 
MAR06710 
MAR06720 
MAR06730 
MAR06740 
MAR06750 
MAR06760 
MAR06770 
MAR06780 
MAR06790 
MAR06800 
MAR06810 
MAR06820 
MAR06830 
MAR06840 
MAR06850 
MAR06860 
MAR06870 
MAR06880 
MAR06890 
MAR06900 
MAR06910 
MAR06920 
MAR06930 
MAR06940 
MAR06950 
MAR06960- 
MAR06970 
MAR06980 
MAR06990 
MAR07000 
MAR07010 
MAR07020 
MAR07030 
MAR07040 
MAR07050 
MAR07060 
MAR07070 
MAR07080 
MAR07090 
MAR07100 
MAR07110 
MARO7120 
MARO7130 
MARO7140 
MARO7150 
MARO7160 
MARO7170 
MARO7180 
MARO? 190 
MARO7200 


Ar 


CC=CINV(IM,IL,IG,T+IT) 
IF(CC .NE.0) THEN 
AR=Y(IM,IL,IG,T+IT)/CC 
ELSE 
AR=0 
ENDIF 
IF(D(IM,IL,IG) .NE. O .AND. PMLE(IM,IL,IG).NE.O 
.AND. PMLE(IM,IL,IG).NE. 1) THEN 


KP=KP+1 
RRML(1T)=RRML(1T)+(CC*(AR-PMLE(1M,1L,1G))**2) 
/(PMLEC IM, IL,1G)*(1) = РИБО СИВИ 
ENDIF 


PRODUCE CHI-SQUARE MEASURE OF EFFECTIVENESS 


RRTSCACIT)= RRTSCA(IT)+CHISQR(DCIM,IL,IG),TSCACIM,IL,IG), 
CC,AR) 
КЕТЅАМ(ІТ)= RRTSAM(IT)+CHISQR(DCIM,IL,IG),TSCAM(IM,IL,IG), 
CC,AR) 
(IT)*CHISQR(D(IM,IL,IG),PJ (IM,IL,IG), 
CC,AR 


RRJ  (IT)2 RRJ 


RRJM (IT)= RRJM (IT)*CHISQR(D(IM,IL,IG),PJM (IM,IL,IG), 
RRJ1 (IT)= RRJ1 С о к И (IM ΤΙ το). 
RRJ1M (IT)= RRJIM CTE és 1D 10 DRE 
RRPEB (IT)= RRPEB (IT) +CHISOR(DCIM, EL, 16) ,PEB (IM,IL,IG), 
RRPTEB (IT)= RRPTEB И ий (IM ID IGA 
RRPEBM(IT)= RRPEBH( IT) +CHISOR(DCIM, IL, 1G) ,PEBM(IM, 1b, 16), 
RRPTBM( IT)= RRPTBMCIT)#CHISQR(DC İM, IL, IG) ,PTEBMCIM, IL, 16), 


RRJLT (IT)2 RRJLT (IT)*CHISQR(D(CIM, IL, IG) ,PJLT( IM, IL, IG), 
CC,AR) 


PRODUCE MEAN ABSOLUTE DEVIATION MEASURE OF EFFECTIVENESS 


IF(DCIM,IL,IG) . МЕ. O .AND. PMLEECIM ΕΡΕ 
AND. PMLE(IM IL IG) ს.) TE 
CALL UOM( CC, AR, PHLE(IM IL IG) MOPMLCECITIS 
MUPMLE(IT), KVOPMLCIT), KVUPMLCITOSSRVBEMDECOIIND 
ENDIF 


IF(CC. NE. 0. ) THEN 
CALL UOM( CC, AR, TSCA(IM,IL,IG), MOTSCA(IT), 
MUTSCAC IT), KVOTSACIT), KVUTSAC IT), KVDTSA(IT)) 
CALL UOM( CC, AR, TSCAM(IM,IL,IG), MOTSAM(IT), 
MUTSAM(IT), KVOTSM(IT), KVUTSM(IT), KVDTSM(IT)) 
CALL UOM( CC, AR, PJ(IM,IL,IG), MOPJ(IT), 
MUPJ(IT), KVOPJ(IT), KVUPJO II აას ОЕ 


CC, AR, PJM(IM,IL,IG), MOPJM(IT), 


CALL UOM( 
* MUPJM(IT), KVOPJM(IT), KVUPJMCIDOSMBUPPIMCOITOM 


CALL UOM( CC, АК, РЈІ(СТМ, ТЕ, 16). МОРО 


T 
aP 
`. 


G4 


MARO7210 
MARO 7220 
MARO7230 
MARO7240 
MARO7250 
MARO7260 
MARO7270 
MARO7280 
MARO7290 
MARO7300 
MARO7310 
MARO7320 
MARO7330 
MARO7340 
MARO 7350 
MARO 7 360 
МАК07370 
MARO7380 
MARO7390 
MARO7400 
MARO7410 
MARO7420 
MARO7430 
MARO7440 
MARO7450 
МАКО 7460 
MARO7470 
MARO7480 
MARO7490 
MARO7500 
MARO7510 
МАК07520 
МАК07530 
МАК07540 
МАК07550 
МАК07560 
MARO7570 
МАК07580 
МАК07590 
МАК07600 
МАК07610 
MARO7620 
MARO7630 
MARO7640 
MARO7650 
MARO7660 
ΜΑΚΟ7670 
MARO7680 
MARO7690 
MARO7700 
MARO7710 
МАК07720 
МАК07730 
МАК07740 
МАК07750 
МАК07760 


270 
260 
250 


240 


CALL UOM( 


CALL UOM( 


CALL UOM( 


CALL UOM( 


CALL UOM( 


ENDIF 


CONTINUE 
CONTINUE 
CONTINUE 


AUN ο απ КИСЕТІСІТУ, КҮШР71(171)) 


ხა გა. PIIMCIM TE, de), HOPJIMCITO; 
πριμ τω. IPSOPIMO PP), KVUP158C 17), KVDPYMCIT)) 


CC, 51, PEB (IM IL IG), MOPEBCIT), 
ყა ა I Ορ. S KVUPEBOIT), KVDPEB(IT)) 


ο PTEBCIMI IL IG) Πο μας το, 
HUPITEBODD). KVOPTBC IT) KVUEEBCIT).. KVDPTBCIT)) 


სა ΡΕΡΗΓΤΗ, ΙΤ. 15}. MOPED II). 
MUPEBM( IT), KVOPBM(IT), KVUPBM(IT), KVDPBM(IT)) 


СС, AR, PIEBMCIM, IL, IG), MOPTBMCIT), 
MUPTBM(IT), КУОВТМ( ІТ), KVUBTM(IT), KVDBTM(IT)) 


mee Kr .EQ. 0) THEN 


FF = 0 
ELSE 


FF-FLOAT(KV(IT))/FLOAT(KP) 


ENDIF 


RRML( IT) ZRRML( IT)*FF 


CONTINUE 


COMPUTE TIME AVERAGE MAD 
MEASURE OF EFFECTIVENESS 


DO 300 IT=1, NVYR 
MUPMLE(IT) = MUPMLE(IT)/KV( IT) 
MUTSCA( IT) = MUTSCA(IT)/KV(IT) 
MUTSAM( IT) = MUTSAM(IT)/KVCIT) 
МОРЈ (ІТ) = МОРЈ (ІТ) /КҮ(ІТ) 
MUPJM (IT) * MUPJM(IT)/KV( IT) 
MUPJ1 (IT) = MUPJ1(IT)/KVCIT) 
MUPJIM(IT) = MUPJIM( IT) /KV(IT) 
МОРЕВ (ІТ) = МОРЕВ(ІТ) /КУ(ІТ) 
MUPTEB (IT) » MUPTEB(IT)/KV(IT) 
MUPEBM( IT) = MUPEBM( IT) /KV(IT) 
MUPTBM( IT) = MUPTBM( IT) /KV(IT) 
MOPMLE(IT) = MOPMLE( IT) /KV(IT) 
MOTSCA(IT) = MOTSCAC IT) /KV( IT) 
MOTSAM(IT) = MOTSAM( IT) /KV( IT) 
MOPJ (IT) = MOPJ (IT)/KVCIT) 
МОРЈМ (ІТ) = МОРЈМ(ІТ) /КУ(ІТ) 
МОРЈ1 (ІТ) = МОРЈ1(ІТ) /КУ(ІТ) 
МОРЈІМ(ІТ) = МОРЈІМ(ІТ)/КҮ(ІТ) 
MOPEB (IT) = MOPEB(IT)/KV(IT) 
MOPTEB (IT) = MOPTEB( IT)/KV(IT) 
MOPEBM(IT) = MOPEBM( IT)/KV(IT) . 
MOPTBM(IT) = MOPTBM(IT)/KV(IT) 
MDPMLE( IT) = MUPMLE( IT) + MOPMLE( IT) 
MDTSCA(IT) = MUTSCA(IT) + MOTSCA(IT) 
MDTSAM(IT) = MUTSAM(IT) + MOTSAM(IT) 


VY 
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MARO7770 
MARO7780 
MARO7790 
MARO7800 


` MARO7810 


MARO7820 
MARO7830 
MARO7840 
MARO7850 
MARO7860 
MARO7870 
МАК07880 
МАК07890 
МАК07900 
МАК07910 
MARO7920 
MARO7930 
MARO7940 
MARO7950 
МАК07960 
МАК07970 
МАК07980 
МАК07990 
MARO8000 
МАК08010 
МАК08020 
MARO8030 
MARO8040 
MARO8050 
МАК08060 
МАК08070 
МАК08080 
МАК08090 
МАК08100 
МАК08110 
MARO8120 
MARO8130 
МАК08140 
МАКО8150 
МАКО8160 
MARO8170 
МАК08180 
МАК08190 
MARO8200 
МАК08210 
МАКО8220 
МАКО8230 
МАК08240 
MARO8250 
МАК08260 
МАК08270 
МАК08280 
МАК08290 
МАК08300 
МАК08310 


MARO8320 - 


300 


ee I A er 


XX 


are TI TFT TFA I kkkt kdekdo terce 


XX 


906 


907 


С ххх 
С Far 


C www 
E A 


С XXX 


908 


Ate 


MDPJ (IT) 
MDPJM (IT) 
MDPJ1 (IT) 
MDPJ1M(IT) 
MDPEB (IT) 


ΠΕ ΠΠ 
Μορ P 
MUPJ1 (IT D; 
MUPJIMCI DET 
МОРЕВ ο. 


MDPTEB (IT) = MUPTIEB (MI) 


HUPEDB TC TD) 
MDETBMCIT) 
CONTINUE 


МОРЕВМ(ІТ) + 
MUPTBM( IT) + 


MOPS (CI 
MOPJM (IT) 
ΡΠ 
MOPJ1M(IT) 
ესს CID 
- МОРТЕВ КЕ 
MOPEBM(IT) 
MOPTBMCIT) 


PRINT SECTION FOR ALL OUTPUT. 


PRINT GROUP CHARACTERISTICS 


WRITE(6,'(1X,''TEN OR LESS BAYES VARIANCE ITERATIONS'')') 


WRITE(6,'(1X,I2,'' BAYES VARIANCE ITERATION(S)'')') IC 
WRITE(6, (1X, XBB Бб BB 


WRITE(6,905) T,KK,SSB,SSE,SHRX,SHR,ZB,ZTB,PZTB,A 
GROUP CHARACTERISTICS , /,' T=',I6, 
SSBs',F12.4, /,' 
SHRS',F12.4, /,' 


FORMAT( 
m 
/» 


la РПІ oa, 
PRINT RESULTS FOR TRANSFORM SPACE 


WRITEC6 >) 


I=1, 


FORMAT( 14X, 10112) 
WRITE(6,907) 'RTSCA ', (RTSCACIT), IT=1,NVYR) 


FORMAT( 
WRITE( 6,907) 
WRITE( 6,907) 
WRITE( 6,907) 
WRITE( 6,907) 
WRITE( 6,907) 
WRITE( 6,907) 
WRITE( 6,907) 
WRITE( 6,907) 
WRITE( 6,907) 


WRITE(6,*) 

WRITE(6,*) ' 
WRITE(6,907) 
WRITE( 6,907) 
WRITE( 6,907) 
WRITE(6, 907) 
WRITE( 6,907) 
WRITE( 6,907) 
WRITE( 6,907) 
WRITE( 6,907) 
WRITE( 6,907) 
WRITE(6,907) 
WRITE( 6,907) 
WRITE(6,907) 
WRITE( 6,908) 


'RSL 
"ЕТ. 


Αι. ορ ДЫ 
"RJ 
'RJ1 
' RXEB 
' RXTEB 
ΕΤ. 
"ЕМІ, | 
! 


pac 
SACK 


(RIJ 


A=’ ხების“ 


' RESULTS FOR TRANSFORM SPACE' 
WRITE COPOS) MANO 


NVYR) 


(IT), IT=1,NVYR) 
(IT), IT=1,NVYR) 


Т (RXEB (IT), IT=1,NVYR) 

, (RXTEB (IT) το αμ 
, (RJLT (IT) ОА) 
, (RML (II) τα 
, (RSL (IT), IT=1,NVYR) 


(IT), IT=1,NVYR) 


'RTJLT ', (RTJLICIS), IT=1,NVYR) 
PRINT RESULTS FOR ORIGINAL SPACE 


RESULTS FOR ORIGINAL SPACE' 


"RRML 


"» (КЕМІ, (ТТ) Е, yp) 


‚AT 
/, ΚΞ 216, 


SSE= ,F12.4, /,' SHRX=',F12.&, 


zB=' ,F12.4, /,' ZTB= ,F12.4, 


AT=" ,F12. 4) 


'RRTSCA', (RRTSCA(IT), IT21,NVYR) 
'RRTSAM', (RRTSAM(IT), IT=1,NVYR) 


ЕЕЗ 
"ВЕЈМ 1 
"RRJ1 ' 
ORR Mes 
'RRPEB ' 
'RRPTEB 
'RRPEBM', 


( RRJ 


=> 7 ә ө м 


(ТТ), ХЕЕЕ) 


(RRJM (IT), IT=1,NVYR) 
(RRJ1 (IT), IT=1,NVYR) 
(RRJ1M (IT), IT=1,NVYR) 
(RRPEB(IT), IT=1,NVYR) 
, (RRPTEB(IT), IT=1,NVYR) 
(RRPEBM(IT), IT=1,NVYR) 


'RRPTBM', (RRPTBMCIT), IT=1,NVYR) 
'RRJLT ', (КЕЛІСІП), ПІ! 


"RV 12: 


(LV 


FORMAT(' ',414,10112) 


WRITE(6,'(9X,''K 


МЕІТЕ(6,%) 


PRINT MAD CALCULATIONS 


CIT), ties) 


''.I12) κΕ 


К; 
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МАКО8330 
MARO8340 
MARO8350 
MARO8360 
МАК08370 
МАК08380 
МАК08390 
МАК08400 
МАК08410 
MARO8420 
MARO8430 
МАК08440 
МАК08450 
МАК08460 
MARO8470 
MARO8480 
МАК08490 
МАК08500 
МАК08510 
МАК08520 
МАК08530 
МАК08540 
МАК08550 
МАК08560 
МАК08570 
МАК08580 
МАК08590 
МАК08600 
МАКО8610 
МАКО8620 
МАК08630 
MARO8640 
МАК08650 
MARO8660 
МАК08670 
МАК08680 
МАК08690 
МАК08700 
МАК08710 
МАК08720 
MARO8730 
МАК08740 
МАК08750 
МАК08760 
МАК08770 
МАК08780 
MARO8790 
МАК08800 
MARO8810 
МАКО8820 
MARO8830 
МАК08840 
МАК08850 
МАК08860 
МАК08870 
МАК08880 


1005 
1006 
1007 


1015 
1016 


ЕЕ СБ (2X РМ’ ')') 
CALL WRMAD(IVYR,NVYR,MUPMLE,KVUPML , MOPMLE , KVOPML , MDPMLE , KVDPML) 
WOSITE(6,' (2X, ' 'TSCA'- )') 
CALL WRMAD(IVYR,NVYR,MUTSCA,KVUTSA ,MOTSCA ,KVOTSA ,MDTSCA,KVDTSA) 
WRITE(6,'(2X,''TSCAM'')') 
CALL WRMAD(IVYR,NVYR,MUTSAM,KVUTSM,MOTSAM , KVOTSM , MDTSAM , KVDTSM) 
WRITE(6,' (2X, ' 'PJ'') ^) 
BALL WRMAD(IVYR,NVYR,MUPJ,KVUPJ,MOPJ,KVOPJ,MDPJ,KVDPJ) 
WRITE(6,'(2X,''PJM'')') 
CALL WRMAD(IVYR,NVYR,MUPJM, KVUPJM , MOPJM , KVOPJM , MDPJM ,KVDPJM) 
WRITE(6,'(2X, pg rt)" ) 
BALL WRMAD(IVYR, NVYR MUPJ1,KVUPJ1,MOPJ1,KVOPJ1,MDPJ1,KVDPJ1) 
WRITE(6,'(2X, 'bJ1M' ) ) 
CALL WRMAD(IVYR,NVYR,MUPJ1M,KVUP1M,MOPJ1M,KVOP1M,MDPJ1M ,KVDP1M) 
WRITE(6,'(2X, ''РЕВ'')' ) 
CALL WRMAD(IVYR,NVYR,MUPEB,KVUPEB , MOPEB, KVOPEB,MDPEB , KVDPEB) 
WRITE(6, (2X, ' ' PEBM' у^) 
CALL WRMAD(IVYR,NVYR,MUPEBM,KVUPBM,MOPEBM, KVOPBM,MDPEBM, KVDPBM) 
WRITE(6,'(2X,''PTEB'')') 
CALL WRMAD( IVYR,NVYR ,MUPTEB , KVUPTB , MOPTEB , KVOPTB , MDPTEB , KVDPTB) 
WRITE(6,'(2X,''PTEBM'')') 
CALL WRMAD(IVYR,NVYR ,MUPTBM,KVUBTM,MOPTBM,KVOBTM,MDPTBM,KVDBTM) 
WRITE(6,*) ! | 
PRINT THE CENTRAL INVENTORY MATRIX 
WRITE(6,*) 'CINV MATRIX' 
DO 1006 IT=1,NYR 
WRITE(6,*) 'YEAR',IT 
DO 1005 IL=1,NLOS 
WRITE(6,1007)(CINVC(IM,IL,1,IT),IM=1,NMOS) 
CONTINUE 
CONTINUE 
FORMAT(15(1X,F7.2)) 
WRITE(6,*) 
PRINT THE ATTRITION MATRIX 
WRITE(6,*) 'Y MATRIX' 
DO 1016 IT=1,NYR 
WRITE(6,”*) 'YEAR',IT 
DO 1015 IL=1,NLOS 
WRITE(6,1007)(YCIM,IL,1,IT) , IM=1,NMOS) 
CONTINUE 
CONTINUE 
WRITE(6,*) Ж 
PRINT ALL ARRAYS 
WRITE(6,'(2X,' "ARRAY TSCA'')') 
CALL ARRAY(MXMOS,MXLOS,MXGRD,NMOS ,NLOS , NGRD,TSCA) 
WRITE(6,'(2X,' ARRAY TSCAM'')!) 
CALL ARRAY(MXMOS,MXLOS ,MXGRD, NMOS , NLOS , NGRD , TSCAM) 
WRITE(6,'(2X,''ARRAY SHREB'')') 
CALL ARRAY(MXMOS ,MXLOS ,MXGRD ,NMOS ,NLOS ,NGRD , SHREB) 
WRITE(6,'(2X,''ARRAY SHRTEB'')') 
GALL ARRAY( MXMOS ,MXLOS , MXGRD , NMOS , NLOS ,NGRD , SHRTEB) 
WRITE(6,'(2X,''ARRAY PMLE'')') 
BALL ARRAY(MXMOS , MXLOS, MXGRD, NMOS , NLOS ,NGRD , PMLE) 
WRITE(6,'(2X, UT ARRAY PJ Sr) 
CALL ARRAY(MXMOS ,MXLOS ,MXGRD ,NMOS ,NLOS ,NGRD,PJ) 
WRITE(6,'(2X,' "ARRAY PJM'')') 


Ta > 
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МАК08890 
MARO8900 
МАК08910 
МАК08920 


` HAR08930 


МАК08940 
МАК08950 
МАК08960 
MARO89 70 
MARO89 80 
МАК08990 
MAR09000 
MAR09010 
MAR09020 
MAR09030 
МАК09040 
МАК09050 
MARO9060 
MARO9070 
MARO9080 
МАК09090 
МАК09100 
МАК09110 
MARO9 120 
МАК09130 
МАК09140 
МАК09150 
MARO9160 
MARO9 170 
МАК09180 
МАК09190 
МАК09200 
МАК09210 
МАК09220 
MARO9230 
МАК09240 
МАК09250 
MARO9260 
МАК09270 
MARO9280 
MARO9290 
МАК09300 
МАК09310 
МАК09320 
MARO9330 
MARO9340 
MARO9350 
MARO9360 
MARO9370 
МАК09380 
MARO9390 
MARO9400 
MARO9410 
MARO9420 
ΜΑΚΟ9 430 
МАК09440 


* 


τος ον ος οσο συ ο κ οσο ης ο ο ο ο ος νο κος κ ολο ο ο ολο νο ο ο ος ο ο κο ο ο κε 


νο ο ο ο ο ο νο ο ο ο ο ο ο ο ο ο ο οσο ο ο ο ο ο ο οσο οσο ο οσο ο ος 


RT TEE TE TC HET TE TC ET TE TC TH TC TH TR RA HT A TH FR TH AR KT HT A RAR KT HT TH AT 


ποσο ο XC XV 2 2C1CXCX XC ededefe vede Jefe CC XC XX XX XX X XC X9C 22920 X 2C X XX XX 2 2ლX X CC 2CXXXCXX CICოXXCICXCXCXXXX 


CALL ARRAY(MXMOS ,MXLOS , MXGRD , NMOS , NLOS , NGRD , PJM) 
WRITE(6,'(2X,''ARRAY PJ1'')') 
CALL ARRAY(MXMOS ,MXLOS , MXGRD , NMOS , NLOS , NGRD , PJ1) 
КЕТТЕ ОБЕ Еа ВИ, 
CALL ARRAY(MXMOS,MXLOS , MXGRD , NMOS , NLOS ,NGRD , PJ1M) 
WRITE(6, (2X,''ARRAY PEB'')') 
CALL ARRAY(MXMOS ,MXLOS ,MXGRD , NMOS ,NLOS , NGRD, PEB) 


WRITE(6, 


' (2X, ' ' ARRAY PTEB'')' 


CALL ARRAY(MXMOS,MXLOS, MXGRD, NMOS, NLOS , NGRD, PTEB) 
WRITE(6,' (2X, ' ' ARRAY PEBM' ')') 
CALL ARRAY(MXMOS ‚MXLOS ‚MXGRD,NMOS ,NLOS ,NGRD, PEBM) 


WRITE(6, 
CALL ARRAY(MXMOS,MXLOS , MXGRD, NMOS , NLOS , NGRD, PTEBM) 


'(2X,'' ARRAY PTEBM'')') 


WRITE(6,'(2X,''ARRAY V'')') 
CALL ARRAY(MXMOS ,MXLOS ,MXGRD ,NMOS ,NLOS ,NGRD,V) 
WRITE(6,'(2X,' "ARRAY XB'')') 


CALL ARRAY (MAMOS , MXLOS , MAGRD , NMOS , NLOS , NORD, XB) 


WRITE(6, 


CALL ARRAY (MXMOS , MXLOS ,MXGRD, NMOS ,NLOS , NGRD, VT) 


"(2X J ARRAY YT OSD 


WRITE(6,'(2X, '' ARRAY XTB'')! ) 
CALL ARRAT(MXMOS, MXLOS , MXGRD , NMOS , NLOS , NGRD , XTB) 


WRITE(6, 


CALL АККАҮ(МХМО5, MXLOS , MXGRD , NMOS , NLOS, NGRD რეი 


WRITE(6, 


2 "ARRAY πα. 


τοσα ! ! ARRAY XJ1'')! ) 


CALL ARRAY (MXMOS , MXLOS , MXGRD ‚NMOS,NLOS,NGRD,XJ1) 


WRITE(6, 


OR ЕВУ 


CALL ARRAY(MXMOS , MXLOS , MXGRD , NMOS ,NLOS , NGRD, XEB) 
WRITE(6, (2X, ' ' ARRAY XTEB'')') 
CALL ARRAY(MXMOS,MXLOS,MXGRD,NMOS ,NLOS ,NGRD,XTEB) 


END 


SUBROUTINE WRMAD(IYR,NYR,MU,KU,MO,KO,MD,KD) 


INTEGER IYR(NYR),KU(CNYR),RO(NYR) , KD(NYR) 


REAL 


MUCNYR) ‚MO(NYR) ,MD(NYR) 


LU CH (Же ' IDEVIATION'' ,T26,12,T42,12,T59 ШШ; 


С (9X, "UNDERAGE! ' το ds 
t 


ж 


ж 
WRITE(6, 
ж 


ж 


RETURN 
END 


( ТҮК(ІМ),ІМ-1,МҮК) 


С ιο ον ‚T55, F9. 6 ¢ 


CES 


">. 


,MUC1), KUC), М0(2), KU( 2), DR KU(3). 
IRITEC6,. (9X, ' 'OVERAGE' ' , 122. pgus (: 


(''.I2,'') T55 55566 SS 


Lx 


‚Mocl), KO(1), μο( 2) .K0(2), M0(3 NE 


' (9X, τ, T22, L9. 6. 5 τς 212, 


(12,48)! 755 О 


3125 


MD(1), KD(1), MD(2), KD( 2), MD(3), KD(3) 
WRITE(6,'(/)') 


SUBROUTINE ARRAY(MM,ML,MG,NM,NL,NG,OBJ1) 


REAL OBJ1(MM,ML,MG) 
DO 10 IG1=1,NG 


` 
4." 
სვ ა 
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t 


δι 


9 
ων 


pi 


Т 
'1 


' ,T38,F9. 6, 
')') 
‚T38,F9.6, 
)') 
38,9. 6, 
) ) 


MARO9450 
МАК09460 
МАК09470 
МАК09480 


' MARO9490 


МАК09500 
МАК09510 
МАК09520 
МАК09530 
МАК09540 
МАК09550 
МАК09560 
MARO9370 
МАК09580 
МАК09590 
MARO9600 
МАК09610 
МАК09620 
МАК09630 
MARO9640 
MARO9650 
MARO9660 
MARO9670 
MARO9680 
MARO9690 
MARO9700 
MARO9710 
MARO9720 
MARO9730 
MARO9740 
MARO9750 
MARO9760 
MARO9770 
MARO9780 
MARO9790 
MARO9800 
MARO9810 
MARO9820 
MARO9830 
MARO9840 
МАК09850 
МАК09860 
МАК09870 
МАК09880 
МАК09890 
МАК09900 
МАК09910 
МАК09920 
МАК09930 
MARO9940 
MARO9950 
MARO9960 
MARO9970 
MARO9980 


MARO9990 · 


DO 10 IL1-1,NL 
10 WRITE(6, (2X,15(F7.3,1X))') (OBJ1(IM1,IL1,1G1), IM1=1,NM) 
WRITE(6, (/) ) 
RETURN 


УУСУ AT A I A A A RA AT ET XXX ee 


BUBROUTINE FTI(CINV,Y, X) 
HAHA Me HAHA HH AAA ICH NICK AICTE HC TIEN FEHR IERIE HIRI TAIT HE HNC HE ETE He TEI He HF EIEN TEN FE 
REAL CINV,X,Y 
TEMP=Y/( 1+CINV) 
ШЕЗРІ-(1. ҒҮ)/( 1“СІКМУ) 
8 P= SORT 5-0 I MV) O A CASINC-1.+2. TEMP) + 
ASIN( <1. +2. *TEMP1) ) 
== 
END 
IHR HT HT HH HIT RICH TH IHR ενος IH IRIN PTR RAIA IR ITA RICH IA IAT ЙС 


SUBROUTINE UOM(VINV,ACT,OBJ,OVER,UNDER,KO,KU,KD) 
Tann N 
REAL OBJ 
ШИСГАСТ -СТ. OBJ) THEN 
UNDER = UNDER+VINV*( ACT-OBJ) 
KU=KU+1 
ELSE 
OVER = OVER +VINV*(OBJ-ACT) 
KO=KO+1 
ENDIF 
KD=KD+1 
RETURN 
END 
TI HII IERIE dede HTT HER TIEN THEA TET რ X дс е Ас с АССА сс С с СТА 


REAL FUNCTION RHO(XB,XBB,SSE,SSB,DEE,KK, T, NN) 
ო_____-_-.+33–--_--_--- ick 

INTEGER T 

REAL*8 SIGSQ, V, U 

SIGSQ=SSE/(NN+2) 

V=SSB/T 

U=KK*(KK-3)*(XB-XBB)**2/(V*(KK-1)) 

RHO=DMIN1(1.DO, DEE/DSQRT(U)) 

END 
ο ου ο ος οι οσο A Ida da ὑς ὑς ος EEE ὑς οσο ο οσο οσους ποσο οσο ο οκ 


REAL FUNCTION MILLER( X, AVGINV ) 
deteiefeiefeiedefetefedfedeiefedededededetse!edeyeyedededededejedefe ife djedjedededefejedesefefejedefedjesejek Jededjededesedede Jet dee fe fedee 


= X/SQRT(AVGINV+. 5)+1. 570796 


TEMP1 = ASIN(SQRT(1. /CAVGINV+1))) 
CIN ПТ. ТЕМР1) THEN 
MILLER = 0.0 
GOTO 1 
ENDIF 
MIM .GT. 3.14159-TEMP1) THEN 
MILLER = 1.0 
GOTO 1 
ENDIF 
ος) . EQ. 0.0) THEN 
MILLER = 0.5 
GOTO 1 
ENDIF 


n 
4 
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MAR10000 
МАК10010 
МАК10020 
МАК10030 


` MAR10040 


МАК10050 
MAR10060 
MAR10070 
MAR10080 
МАК10090 
МАК10100 
МАК10110 
МАК10120 
MAR10130 
MAR10140 
МАК10150 
МАК10160 
МАК10170 
МАК10180 
МАК10190 
МАК10200 
MAR10210 
MAR10220 
MAR10230 
MAR10240 
MAR10250 
MAR10260 
МАК10270 
МАК10280 
МАК10290 
МАК10300 
МАК10310 
МАК10320 
МАК10330 
MAR10340 
MAR10350 
MAR10360 
МАК10370 
МАК10380 
МАК10390 
MAR10400 
MAR10410 
MAR10420 
MAR10430 
MAR10440 
MAR10450 
MAR10460 
MAR10470 
MAR10480 
MAR10490 
MAR10500 
MAR10510 
MAR10520 
MAR10530 
MAR10540 
MAR10550 


IE (COS(TM) .GT. 0.0) = 1 
IF (COS(TM) . LT. 0.0) ვი. l 
TEMP2=1-(SINCTM)+CSINCTM) -(1. /SINCTM) ) ) /AVGINV)**2 
ТЕ (TEMP2. LT. 0.0) TER 2080 
MILLER=. 5*( 1. -SIGNUM*TEMP2**. 5) 
N 
И. 


REAL FUNCTION SCINV(X,A) 
Tee 


XX SCALE INVERSION (TRANSFORM ==> ORIGINAL SCALE) 
R=X/SQRT(A+. 5) 
ТЕ(В .ТТ. -1. 570796) THEN 
SCINV=0. 
ELSE IF(R .GT. 1.570796) THEN 
SCINV=1. 
ELSE 
SCINV=. 5*( 1. +SIN(R)) 
ENDIF 


λος οσους ος ο Weve seve Fete Fes TE Te FT АО 


REAL FUNCTION VAR( ZZ ) 
3e ese ese esee Sede desee ede dese yedede desee de see еее λος κος ο ΧΧΧ οσο κος ο ο κος 
DATA AA/1.6835/, B1/-.8934/, B2/.8991/ 
IL C 2250) EN 
VAR = 1. 
ELSE 
IL ( 772 στ |, ეცი ШЕН 
VAR = .05 
ELSE 
VAR = AA*(ZZ**B1)*(22-1. )**B2 
ENDIF 
ENDIF 
END 
ee ee ee ee ee de reed Te ہاو لو او‎ Jede te fe de tele Tele Je Je Ve Sec ke Ich Tele le tete rettete te Fe tesede 


REAL FUNCTION CHISQR(ID,OBJ,INV,ACT) 
ολο ο ο οσο ο ο Jose ο ος ος ο πμ ποπ πμ κκ πμ πι. - 
REAL OBJ,INV 
IF (ID.NE.0 . AND. ОВЈ. МЕ. О. . AND TOBTSNE SIRT) 
* CHISQR = (INV*(ACT-0BJ)**2)/(OBJ*(1. -0BJ)) 
END 


seqe οσο) ο ο ολ ολο κος ος λος ος οσο ος κος κο ο κο ο ΧΧΧ ΝΟ eFC Fe 


SUBROUTINE READER(MOS,LOS,GRADE,YR,INV,Y, NMOS,SMOS, 


*  NGRD,SGRD, 51051,51052, 511,512, ШІП ІСТЕУШІ 
HH HEHE LUTHER RTL HT TRTRLU TUT TR RT KUH TA HH A AT AH AT HT 


ο READ RECORD AND ACCUMULATE LOSSES 
INTEGER MOS, LOS, GRADE, YR, INV, LOSS(8) 
REAL Y 


ІМТЕСЕК 5Т1, ST2 
INTEGER SLOS1, SLOS2 
INTEGER SMOS(40), NMOS 
INTEGER SGRD(10), NGRD 


ICNT=0 : А 
1 READ(1,100,END=999) YR,MOS,GRADE,LOS,INV,(LOSS(I),I=1,8) 
100 ЕОКМАТ( 412,915) 
XX CHECK IF RECORD MEETS SELECTION CRITERIA. OTHERWISE REJECT. 
ІМ-0 


59 
x... 
è 
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т 


МАК10560 
МАК10570 
МАК10580 
МАК10590 
МАК 10600 
МАК10610 
МАК10620 
МАР.10630 
МАК 10640 
МАК10650 
МАК 10660 
МАК10670 
МАК 10680 
МАК10690 
МАК10700 
МАК 10710 
МАК10720 
МАК10730 
MAR10740 
MAR10750 
MAR10760 
MAR10770 
MAR10780 
MAR10790 
MAR10800 
MAR10810 
MAR10820 
MAR10830 
MAR10840 
МАК10850 
МАК10860 
МАК10870 
МАК10880 
МАК10890 
МАК10900 
МАК10910 
МАК10920 
МАК10930 
МАК10940 
МАК10950 
МАК10960 
МАК10970 
МАК10980 
МАК10990 
МАК11000 
МАК11010 
МАК11020 
МАК11030 
MAR11040 
МАК11050 
МАК11060 
МАК11070 
МАК11080 
МАК11090 
МАК11100 


MAR11110 · 


10 


20 


30 
40 


999 


DO 10 I= 


TC MOS ΕΠ SHOS(1)) THEN 


END IF 
CONTINUE 
ES TO I 


CONTINUE 
DO 30 I- 


ENSGGRADE .EQ- SGRD(I)) THEN 


END IF 
CONTINUE 
GO TO 1 
CONTINUE 


M LCOS LT. SEOSI .OR. LOS .GT. SLOS2) GO TO 1 


TL 95-5 


YEARS OVER ST2 ARE USED FOR VALIDATION 
ШЕГҮК .LT. ST1) GO TO 1 


IT-YR-ST 


COMPUTE TOTAL LOSS 


Ξ0 
DO 50 I= 
=Y + L 
CONTINUE 


RETURN 
IEOF=1 
END 


1,NMOS 


IM-I 
GO TO 20 


1,NGRD 


IG-I 
GO TO 40 


LOS1+1 


IST 


1,8 
OSS(I) 


x A. 
№ 
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МАК11120 
МАК11130 
MAR11140 
МАК11150 


“МАК11160 


МАК11170 
МАК11180 
МАК11190 
МАК11200 
МАК11210 
МАК11220 
МАК11230 
MAR11240 
MAR11250 
МАК11260 
МАК11270 
МАК11280 
MAR11290 
MAR11300 
МАК11310 
МАК11320 
МАК11330 
MAR11340 
MAR11350 
MAR11360 
МАК11370 
МАК11380 
МАК11390 
МАК11400 
MAR11410 
MAR11420 
MAR11430 


APPENDIX C. SIMULATION FOR CHOICE OF AVERAGE 


PROGRAM MCSIM 
Tee 


tet ete tt eae + + + + + + + + + x + FH + + +* + + + + + 


INVENTORY VALUES 


SIMULATION COMPARISON OF INVENTORY HEAN PERFORMANCE 
PURPOSE: TO IDENTIFY THE BEST MEAN VALUE (ARITHMETIC, 


GEOMETRIC AND HARMONIC) FOR USE IN THE INVERSION OF 
TRANSFORMED DATA TO THE ORIGINAL SPACE. 


LAMBDA - 
TRIALS - 
REPS = 
BASIC - 
ЕТЕ = 
ЕТЕ? = 


7 დ 
MEANS - 


IRSUN = 
TRMULT = 


TRINV - 


INPUT/OUTPUT: 


PROGRAMMER: CAPT C R DICKINSON, USMC 
VARIABLE DESCRIPTION: 


PARAMETER USED IN THE POISSON RANDOM NUMBER 
GENERATOR 

VECTOR OF INVENTORY VALUES GENERATED FROM 
IHE POISSON RANDOM NUMBER GENERATOR 

NUMBER OF INVENTORY VALUES TO BE GENERATED 
BASIC INVERSION OF TRANSFORMED VALUES 
FREEMAN-TUREY EXACT INVERSION WITH 


FREEMAN-TUREY EXACT INVERSION WITH 


TRANSFORMED SCALE VARIABLE 

VECTOR OF ARITHMETIC, GEOMETRIC AND НАКМОМІС 
MEANS 

SUMMATION OF TRIALS FOR USE IN THE 
ARITHMETIC MEAN 

MULTIPLICATION OF TRIALS FOR USE IN THE 
GEOMETRIC MEAN 

SUMMATION OF THE INVERSE OF TRIALS FOR USE 
IN THE HARMONIC MEAN 


INPUT - HARD CODED PARAMETERS 


OUTPUT = SIM OUI A1 (DIST TIT OY 


Uo Gb Cb GB Gb Gb Gt Ct Gto Gto GR GR Gb Gb Gb Gb Gb Gb Ob bo GRoGR o0 o0 oL oX 


+ 


+ 


ж 


Зое УСС С ОСУУ СУУСУ ССС СУС УС УСС УСС УСУ САС ο πω ΚΚ κκ ο ο κ κκ ολ VE TE VEE $ 


e 
oe 


VARIABLE DECLARATION 


INTEGER TRIALS(100),Y(100),REPS(3),PROB,R, 
G COUNT, TRSUM, YTEMP, YSUM 

REAL 2(100),T(100),MEAN(3),BASIC(3),DIFF(3),LAMBDA(5), 
C P(#4),FTE1(3),FTE2C( 3) ოფ, 

DOUBLE PRECISION DPSEED 


CALL EXCMS('FILEDEF 02 DISK REPLIC LISTING Al') 


PI=3. 14159 
DPSEED=889246.DO 


LAM=5 
R=3 
PROB=4 
COUNT=2 


pis 


VARIABLE INITIALIZATION 


1% 


x A 
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MCS00010 
MCS00020 
MCS00030 
MCS00040 
MCS00050 
MCS00060 
MCS00070 
MCS00080 
MCS00090 
MCS00100 
MCS00110 
MCS00120 
MCS00130 
MCS00140 
MCS00150 
MCS00160 
MCS00170 
MCS00180 
MCS00190 
MCS00200 
МС500210 
МС500220 
MCS00230 
MCS00240 
MCS00250 
MCS00260 
MCS00270 
MCS00280 
MCS00290 
MCS00300 
MCS00310 
MCS00320 
MCS00330 
MCS00340 
MCS00350 
MCS00360 
МС500370 
МС500380 
МС500390 
МС500400 
MCS00410 
MCS00420 
MCS00430 
MCS00440 
MCS00450 
МС500460 
MCS00470 
MCS00480 
MCS00490 


XXX 


XXX 


triste 


C 
C 


DATA REPS/20,40,60/ 
DATA LAMBDA/1. ‚2. ‚4. ‚8. ‚16. / 
DATA P/.05,.1,.2,.4/ 


HEADER PRINT 
ΕΙΤΕ (02.'(/T6,' LAMBDA'',T15, 'REPS'',T27, 'P'', 
T34,''METHOD'' ,T46,''ARITHMETIC'' ,T59, 

'"'GEOMETRIC'',T71, 'HARMONIC' !/)') 


LOOP FOR LAMBDA PARAMETER 

DO 500 I=1,LAM 
DO 450 J=1,R 
CALL GGPOS( LAMBDA( 1) ,DPSEED,REPS( J) , TRIALS, IER) 


NSUM=0 
DOTIO I]=1,REPS(J) 


TRIALS(J1)=TRIALS(J1)+1 


NSUM=NSUM+TRIALS( J1) 


LOOP FOR PROBABILITY PARAMETER 
DO 400 K=1,PROB 
COUNT=COUNT+4 
LINE COUNTER FOR PAGE BREAK 
IF (COUNT. GT. 55) THEN 
WRITE (2,'(1H1)') 
ONT 2 
WRITE (02,'(/T6,''LAMBDA'' 
T34,''METHOD'' ,T46, ' 


} 


ШЕР “БЕРӘ А T27. RD 
ARITHMETIC'' ,T59, 


''GEOMETRIC'',T71, 'HARMONIC' '/)'!) 


ENDIF 


TRSUM=0 
TRINV=0 


TRMULT=0. 


ARCSUM=0. 
TSUM=0. 
YSUM=0 


LOOP FOR REPITITION PARAMETER 
DO 20 K1=1,REPS(J) 


CALL GGBN(DPSEED,1,TRIALS(K1),P(K),YTEMP) 


YSUM=YSUM+YTEMP 
TRSUM=TRSUM+TRIALS(K1) 


TRMULT=TRMULT+ALOG( REAL( TRIALS(K1))) 
TRINV=TRINV+1. /REAL(TRIALS(K1)) 


Y(K1)=YTEMP 


Z(K1)=. 5*(TRIALS(K1)+. 5)**. 5*(ASIN( 2. *Y(K1)/ 
(TRIALS(K1)+1. )-1. )+ASINC2.*(Y(K1)+1.)/ 
(TRIALS(K1)+1. )-1. )+PI) 

T(K1)=Z(K1)/( (TRIALS(K1)+. 5)**. 5) 


ARCSUM=ARCSUM+Z(K1) 
TSUM=TSUM+T(K1) 


ZAVG-ARCSUM/REAL(REPS( J) ) 


TAVGZTSUM/REAL(REPS( J) ) 


MEAN VALUE CALCULATION 


E 
1 
„А. 
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MC500500 
MC500510 
MCS00520 


µου 


‘MCS00540 


MCS00550 
MCS00560 
MCS00570 
MCS00580 
MCS00590 
MCS00600 
MCS00610 
MCS00620 
MCS00630 
NCS00640 
MCS00650 
MCS00660 
MCS00670 
МС500680 
MCS00690 
MCS00700 
MCS00710 
MCS00720 
MCS00730 
MCS00740 
MCS00750 
MCS00760 
MCS00770 
MCS00780 
МС500790 
MCS00800 
MCS00810 
MCS00820 
MCS00830 
MCS00840 
MCS00850 
MCS00860 
MCS00870 
MCS00880 
MCS00890 
MCS00900 
MCS00910 
МС500920 
MCS00930 
MCS00940 
MCS00950 
MCS00960 
MCS00970 
MCS00980 
MCS00990 
MCS01000 
МС501010 
MCS01020 
MCS01030 
MCS01040 
MCS01050 


30 


40 


400 
450 
500 


MEAN( 1)=TRSUM/REAL(REPS( J) ) 
MEAN(2)=EXP(TRMULT*(1. /REAL(REPS(JI)))) 
MEAN(3)=REAL(REPS(J))/TRINV 


IF (COS(TAVG).GE.0.0) THEN 
SIGNUM=1. 

ELSE 
SIGNUM=-1. 

ENDIF 


DO 30 K2=1,3 

TAVG1=ZAVG/(MEAN(K2)+.5)%*.5 

IF (COS(TAVG1).GE. 0. 0) THEN 
SIG1=1. 

ELSE 
5161=-1. 

ENDIF 

BASIC AND FTE INVERSE METHOD CALCULATION 

FTEI(K2)=. 5*(1. -(SIGNUM*(1. -(ABSCSINCTAVG)+ 

(SIN(TAVG)-(1. /SIN(TAVG)))/MEAN(K2)))**2. )**.5)) 


FTE2(K2)=. 5*(1. -(SIGI*(1. -(ABS(SIN(TAVG1)+ 
(SIN(TAVG1)-(1. /SINCTAVG1)))/MEAN(K2)))**2. )**, 5)) * 


BASIC(K2)=. 5*( 1. +SINCZAVG/((MEAN(K2)+.5)**., 5) 
ΙΙ 
OUTPUT STATEMENTS 
DO 40 2 
IF (K3. EQ. 1) THEN 
МЕІТЕ (02, (Т6;,Ғ8.1,715,13,125 55.5 785, "ТЕ 
T&6 ,F7.& TS9 F7 6 7 T71 F7 65) 
LAMBDA(I),REPS(J),P(K),(BASIC(IK),IK=1,3) 
ELSE | 
IF (K3. EQ. 2) THEN 
WRITE (02, '(1T35,' FUE] το ҰТТЫ IN 
T71,F2.4) ) (FTE1(IK),IK=1,3) 
ELSE 
ΝΕΙ͂ΤΕ (02, (Τ35. "ЕТЕ? 1546. 87 А ТЗЭ 
T71,F7.4/)') (FTE2(IK),IK=1,3) 
ENDIF 
ENDIF 
CONTINUE 


CONTINUE 
CONTINUE 
CONTINUE 
STOP 
END 


: А, 
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MCS01060 
MCS01070 
MCS01080 
МС501090 
МС501100 
МС501110 
МС501120 
МОЕ О 
МС501140 
МС501150 
МС501160 
МС501170 
МС501180 
MCS01190 
MCS01200 
MCS01210 
MCS01220 
MCS01230 
MCS01240 
MCS01250 
MCS01260 
МС501270 
МС501280 
МС501290 
MCS01300 
MCS01310 
MCS01320 
MCS01330 
MCS01340 
МС501350 
МС501360 
МС501370 
МС501380 
МС501390 
MCS01400 
MCS01410 
MCS01420 
MCS01430 
MCS01440 
MCS01450 
MCS01460 
MCS01470 
МС501480 
MCS01490 
MCS01500 
MCS01510 
MCS01520 
MCS01530 
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